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Introduction 

Since  the  introduction  of  the  micromachining  process,  wherein  mechanical 
structures  are  etched  from  blocks  of  silicon,  a  number  of  very  small  acoustic 
and  vibration  sensors  have  been  built  (for  example,  [1-7]).  Ranging  from 
simple  capacitive  pressure  sensors  [8]  to  accelerometers  that  measure  the 
proof-mass  displacement  by  electron  tunneling  with  active  mechanical  feedback 
[9],  these  sensors  are  attractive  for  many  space-limited  applications. 

However,  the  small  moving  parts  are  especially  susceptible  to  mechanical  noise 
resulting  from  molecular  agitation.  If  the  sensor  Is  Intended  for  low-level 
signals,  this  mechanical-thermal  noise  can  be  the  limiting  noise  component. 

Because  mechanical-thermal  noise  has  not  been  important  for  conventional¬ 
ly  sized  sensors,  the  analysis  procedures  are  often  unfamiliar  and  the 
mechanism  is  often  neglected.  Many  designers  fall  into  the  trap  of  making 
elaborate  analyses  of  more  common  but  less  important  sources  of  noise.  The 
purpose  of  this  paper  is  to  present  some  simple  techniques  for  evaluating 
mechanical-thermal  noise  limits  and  to  review  shot,  1/f,  and  amplifier  noise. 
While  these  techniques  arc  especially  useful  for  designing  microminiature 
sensors,  the  principles  are  generally  applicable  for  any  sensor  Intended  for 
extreme  sensitivity. 

The  basic  types  of  noise  considered  are: 

Thermal  Noise:  The  random  fluctuations  result  from  molecular  vibration. 
This  component  is  called  Johnson  Noise  in  electrical  systems  and 
Brownian  Motion  in  mechanical  systems. 

Shot  Noise:  The  noise  associated  with  random  emission  of  particles  (or 
photons)  that  have  not  yet  reached  thermal  equilibrium  with  their 
surroundings.  Examples  are  emission  of  photons  from  a  laser  and  transit 
of  charge  carriers  across  a  semiconductor  Junction. 
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1/f  Noise:  A  frequently  observed  noise  component  with  power  that  decreases 
roughly  linearly  with  Increasing  frequency  (roughly  equal  power  per 
octave).  As  yet,  there  Is  no  satisfactory  physical  model  for  1/f  noise. 
This  type  of  noise  makes  DC  and  very  low  frequency  measurements  diffi¬ 
cult.  It  Is  also  called  flicker  noise  or  excess  noise. 

Since  both  mechanical,  electrical,  and  thermodynamic  quantities  will  be 
used  In  these  notes,  any  system  of  units  other  than  SI  would  make  calculations 
awkward  and  prone  to  error.  Consequently,  SI  units  will  be  used  throughout: 
temperature  In  kelvln  (K),  mass  In  kilograms  (kg),  and  length  In  meters  (m). 

One  point  of  caution:  these  units  lead  to  a  base  unit  of  pascal  (Pa)  for 
pressure  while  the  standard  reference  for  underwater  acoustics  Is  the  micro- 

pascal  (pPa).  Also,  a  common  unit  for  acceleration  Is  the  "g"  but  the  SI  unit 

2  2 
Is  m/s  ;  one  g  Is  9.8  m/s  . 

All  of  these  noise  quantities  are  distributed  In  frequency  so  they  are 

usually  given  In  terms  of  spectral  density:  some  quantity  (proportional  to 

power)  per  hertz.  Some  examples  of  the  units  are:  voltsVhz,  mVHz.  g^/Hz, 

2 

pPa  /Hz.  These  are  frequently  expressed  In  terms  of  the  linear  quantities  as 
follows:  volts/v'liz,  m/vISz,  g/^Hz,  pPa/i/Hz.  These  linear  units  should  be  used 
with  caution  because,  when  Individual  noise  sources  are  Incoherent,  as  Is  most 
often  the  case,  theli-  powers  add,  not  their  voltages.  To  adjust  for  a  known 
frequency  band,  multiply  the  power-like  quantity  (voltsVHz,  for  example)  by 
the  bandwidth;  or,  multiply  the  linear  quantity  (volts/i/Hz,  for  example)  by 
the  square  root  of  the  bandwidth. 

In  order  for  a  noise  analysis  to  have  practical  value,  the  results  need 
to  be  compared  to  the  desired  sensitivity  of  the  sensor.  In  underwater 
acoustics,  a  reasonable  choice  for  the  desired  minimum  signal  to  be  measured 
is  the  quietest  level  of  ocean  background  noise.  Table  I  gives  an  approxima¬ 
tion  to  these  levels  from  0.1  to  1000  Hz.  The  values  were  obtained  by 
subtracting  10  dB  from  an  average  quiet-ambient  curve  [10].  Notice  that  the 
acoustic  displacements  are  very  small;  the  displacement  at  1000  Hz  Is  about 
equal  to  the  diameter  of  a  proton!  The  pressure  and  the  acceleration  corres- 
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ponding  to  this  SSO-lOdB  ("sea-state  zero  minus  10  dB")  curve  are  shovm  in 
Figure  1  and  tabulated  in  Table  I.  In  practice,  each  sensor  application  must 
be  examined  individually  to  determine  the  minimum  sign?,!  retirement;  in  these 
notes,  however,  the  SSO-lOdB  curve  will  be  used  to  Illustrate  the  analysis. 


I 

M 


Z 

o 

I 


Figure  1.  Acoustic  pressure  anJ  particle  acceleration  correapondlng  to 
Sea  State  Zero  ninut  10  dB. 


Tible  I.  Sea  State  Zero  - 10  dB 


FREQUENCY 

LEVEL 

PRESSURE 

VELOCITY 

DISPLACEMENT 

ACCELERATION 

(f) 

(NL) 

(p) 

(V) 

(*) 

(•) 

Hz 

dB 

uPapiH 

n/a  prH 

m  prH 

nlt^ZvM 

0.1 

122 

12600-)0 

8.4E-07 

l.:,B-06 

53E-07 

1 

92 

400CI 

2.7E-08 

4J£-09 

1.7E4)7 

10 

62 

1260 

8.4E-10 

13E-11 

SJ£4)8 

100 

46 

210 

13E-10 

ZlE-13 

8.4E-08 

1000 

34.5 

53 

3.5E-11 

3.6E15 

Z2E-07 

Notes:  Level  (NL)  ia  in  db//uPa^  2  per  hertz 


"prH"  standa  for  per  root  rertz 
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(The  various  acoustic  amplitudes  are  related  as  followii: 

p  «  10*^“  in  pPa/vliz 

V  «  (p/10®)/pc  (note  conversion  to  Pa; 

pc  “  1  500  000  kg/m^s;) 

X  a  v/«  where  u  *  2nf 

a  a  V(d  ) 

For  audio  microphone  design,  ambient  levels  are  often  expressed  as 
A-weighted  levels  [11]  (where  the  reference  level  is  20  pPa  not  1  liPdi) .  The 
equivalent  bandwidth  for  a  source  of  white  noise  (uniformly  distributed  in 
frequency)  is  roughly  15  000  Hz;  noise  spectral  density  values  can  be 
converted  to  A-welghted  levels  by  multiplying  the  per  hertz  value  by  15  OOOHz 
and  then  dividing  by  400  (20  pPa  squared).  If  the  noise  spectral  density  it 
given  in  dB  with  respect  to  1  pPa^  per  hertz,  the  value  can  be  converted  to 
A-weighted  level  by  adding  16  dB.  Some  typical  background  levels  [11,  12] 
are:  25-30  dB(A)  in  recording  studios,  35-40  dB(A)  in  churches,  and  30-50 
dB(A)  in  quiet  residential  areas 


Thermal  Noise 

The  concept  of  thermal  equilibrium  [13]  is  generally  taught  in  basic 
thermodynamics  courses:  A  collection  of  molecules  reaches  equilibrium  in 
which  each  molecule  has,  on  average,  the  same  amount  of  energy.  If  another 
molecule  is  added  to  the  "bath,"  it  too  reaches  the  same  average  energy  level 
after  a  few  collisions.  What  is  not  often  appreciated  Is  that  this  uniform 
distribution  of  energy  applies  to  macroscopic  objects  also.  A  golf  ball 
placed  in  that  molecular  bath  will,  after  many  collisions,  acquire  an  average 
klneti'.;  energy  equal  to  that  of  any  of  the  molecules.  The  only  difference  is 
that  it  takes  a  substantially  greater  number  of  collisions  to  bring  the  golf 
ball  into  thermal  equilibrium  [14]. 

One  of  the  more  well  known  mechanisms  for  mechanical -thermal  noise  is 
Brownian  Motion  [15].  Here,  the  agitation  of  an  observable  object  is  caused 
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by  mol«cuIar  collisions  from  a  surrounding  gas  or  liquid  and  the  agitation  is 
directly  related  to  the  fluid’s  viscosity.  In  fact,  any  molecular  agitation 
even  through  solid  structures  like  springs  and  supports  can  cause  movement  of 
an  object  and  each  fluctuation  component  is  related  to  a  mechanical  damping. 


Compare  the  lossless  harmonic  oscillator  with  the  damped  harmonic 
oscillator  in  Fig.  2. 


LOSSLESS 


DAMPED 


/////////// 

I 

z 


////////////// 
- ! - 1 - 


m  +  kz  =  0  m  R  —  +  kz  =  0 

dt  dt^  dt 

Flguro  2.  Loiiial«sa  and  damped  simple  harmonic  oaclllatora. 

The  presence  of  damping  in  the  system  on  the  right  suggests  that  any  oscilla¬ 
tion  would  continue  to  decrease  in  amplitude  forever.  Even  the  small,  random 
Jitter  caused  by  molecular  motion  would  decay;  this,  however,  is  the  same 
thing  as  saying  that  the  temperature  (a  measure  of  this  molecular  motion) 
would  continually  decrease.  The  temperature  cannot  drop  below  that  of  the 
surroundings  so  the  model  on  the  right  violates  the  Second  Law  of  Thermo¬ 
dynamics.  For  mechanical  systems,  this  is  not  normally  a  problem  but,  if  a 
sensor  is  being  designed  for  very  small  signals,  the  excitation  of  the 
oscillator  by  molecular  vibration  can  be  significant. 

The  correct  physical  model  for  the  damped  simple  harmonic  oscillator 
Includes  a  force  generator  of  sufficient  amplitude  to  bring  the  system  into 
thermal  equilibrium  rather  than  to  let  it  decay  to  the  equivalent  of  absolute 
zero  temperature.  The  proper  differential  equation  is: 
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m  ^  4  R  ♦  k  z  =  f  (t)  (1) 

j*.  2  n 


where  f  Is  the  fluctuating  force  required  to  maintain  thermal  equilibrium, 
n 

The  presence  of  the  damping  term  in  the  equation  requires  that  the  fluctuating 
force  be  present  as  well.  This  is  a  crucial  relationship. 

Since  noise  sources  add  Incoherently  (power  adds,  not  force),  the  forcing 
power  should  be  linear  with  R.  Also,  the  energy  in  molecular  vibrations  is  a 
linear  function  of  temperature,  so  the  forcing  power  should  also  be  linear 
with  T.  Consequently,  the  force  should  be  proportional  to  the  square  root  of 
R  times  T. 

The  association  between  damping  and  fluctuation  is  expressed  by  the 
Fluctuation-Dissipation  Theorem  [16]:  any  mechanism  for  dissipation  causes 
fluctuations.  If  there  is  a  path  for  energy  to  leave  an  object  or  system, 
then  that  path  also  allows  thermally  excited  molecules  in  the  enviroiunent  to 
transfer  disordered,  fluctuating  energy  to  the  object.  This  "dissipation''  can 
represent  any  mechanism  that  would  cause  the  oscillation  to  decay:  mechanical 
damping,  acoustic  or  electromagnetic  radiation,  transfer  of  power  to  a  load. 

The  basic  relationship  that  governs  thermal  noise  analysis  is  expressed 
by  the  Equipartition  Theorem  [13,  14,  17]:  any  "mode"  of  a  system  in  thermal 
equilibrium  has  an  average  noise  energy  of  T  where  k  is  Boltzmann’s 
constant  (1.38  x  lO"^^  J/kelvln)  and  T  is  the  absolute  temperature  (here  in 
kelvin:  K  »  *C  4-  273) .  The  modes  of  the  system  can  be  identified  by 
writing  out  all  of  the  components  of  the  system  energy.  Any  component  that 
depends  on  the  square  of  a  coordinate  is  a  mode  in  the  above  sense.  Common 

examples  include  kinetic  energy  (-'-mv*),  spring  potential  energy  (^kx^),  energy 

1  2  *  1*2 
stored  in  a  capacitor  (-CV  ),  and  energy  stored  in  an  inductor  (-L1  ). 

Equipartit'iOii  claims  that  the  thermal  energy  is  equally  distributed  among  all 

the  energy  storage  modes  and,  furthermore,  is  equal  to  ^k  T  for  each  mode. 

2  8 
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The  key  to  analyzing  mechanical  systems  in  equilibrium  is  that  simple 

equlpartition  holds  for  every  mode  of  energy  storage  (as  long  as  k  T  »  hf 

"34  ^ 

where  h  is  Plank’s  constant.  6.6  x  10  Js.  and  f  is  frequency).  Each  mode 
(kinetic  or  potential),  whether  corresponding  to  an  observable  motion  of  a 
macroscopic  mass  or  microscopic  motion  of  molecular  motion,  vibration,  or 
rotation,  has  the  same  amount  of  thermal  energy  and  these  modes  are  continual¬ 
ly  exchanging  energy  with  all  of  the  other  modes  in  order  to  maintain  this 
equilibrium. 

If  an  "ordered"  mode  such  as  the  vibration  of  a  macroscopic  mass-spring 
system  is  in  equilibrium  with  a  large  number  of  Independent  microscopic  modes 
(say,  individual  molecular  vibrations  in  the  spring),  then  an  observable 
"dissipation"  results  (IS).  Energy  from  the  orderly  motion  of  the  mass-spring 
system  is  transferred  (at  some  measureable  rate)  to  a  very  large  number  of 
molecular  modes.  In  turn,  the  molecular  vibrational  modes  transfer  energy 
back  to  the  mass-spring  oscillator;  however,  since  the  molecular  vibrations 
are  essentially  independent  of  each  other,  the  return  transfers  are  random  and 
the  probability  of  them  being  in  phase  enough  to  reinforce  the  mass-spring 
motion  is  effectively  zero.  So  the  whole  system  (mass-spring  vibrator  and 
molecules  in  the  spring),  if  no  longer  driven  externally,  eventually  reaches  a 
state  where  the  energy  in  the  mass-spring  motion  equals  the  energy  in  any  of 
the  many  molecular  vibrations. 

The  principle  of  equlpartition  leads  directly  to  the  first  technique  for 
analyzing  mechanical-system  noise  [19].  The  average  displacement  of  the  mass 
in  a  mass-spring  oscillator  is  given  by 


1  1.  ^  2^ 

-  k  <  X  > 

1 

k  T 

2 

i 

B 

Here,  <  >  denotes  an  average  and  it  is  a  broadband  average.  Also,  k  is  the 
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spring  constant  while  k  Is  Boltzmann's  constant.  The  average  velocity  of  the 

B 

mass  is 


B 

A 

< 

10 

V 

1 

k  T 

2 

2 

B 

(3) 


A  free  particle  (an  atom  of  Helium  gas,  for  example)  may  not  have  a 
spring-like  storage  mode  but  its  kinetic  energy  is  still  subject  to  equipar- 
tition.  Since  the  throe  coordinate  directions  are  independent,  there  are 
three  kinetic  modes.  The  x,  y,  and  z  components  each  have  an  average  kinetic 

12  1  3 

energy  of  -m<v  >  =  -k  T  so  the  total  average  energy  of  the  particle  is  -k  T. 

2  X  2  B  2  B 


In  some  analyses,  these  simple  relations  will  suffice  but,  often,  the 
frequency  distribution  of  the  noise  is  required.  In  this  case,  the  noise 
response  of  a  device  can  be  calculated  by  inserting  a  force  generator  at  the 
same  location  as  each  damper  in  the  system.  The  spectral  density  of  the  force 
corresponding  to  each  damper  is  (see  Appendix  A) 


I - 

F  -  v'SiTTO  [N/VHz] 
n  B 


(4) 


This  result  ij  derived  directly  from  equipartltion:  a  system  mode  having  ^k  T 
broadband  energy  is  equivalent  to  the  damper  having  an  associated  force 
generator  with  /4k  TR  spectral  density.  (The  damping  coefficient  R  is 

B 

mechanical  resistance  —  force  per  velocity  —  and  can  be  a  function  of 
frequency.  To  analyze  electrical  systems,  put  a  voltage  generator  of  /4k  tR 

B 

volts//Rz  in  series  with  each  resistor  in  the  circuit.  Here,  R  is  electrical 
resistance  [20,  21].) 


One  of  the  consequences  of  the  Fluctuation-Dissipation  theorem  is  that, 
if  the  dissipation  (loss,  damping,  resistance)  of  a  system  is  measured,  the 
individual  contributors  to  that  dissipation  need  not  be  known;  it  is  the  total 
dissipation  acting  on  a  mass-spring  oscillator  that  determines  the  thermal 
fluctuations  of  that  oscillator.  This  is  why  a  measurement  of  the  Q  of  an 
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oscillatory  system  Is  useful  in  noise  analysis:  the  measured  Q  incorporates 
all  of  the  dissipation  mechanisms  that  are  acting  on  the  system. 

Another  consequence  of  the  F-D  relations  is  that  individual  dampers  can 
be  treated  individually  if  that  is  more  convenient.  This  is  particularly 
useful  in  modeling  devices.  Construct  the  equivalent  circuit  for  the  entire 
device  (include  both  mechanical  and  electrical  parts)  and,  for  each  resis¬ 
tance,  add  a  series  force  (or  voltage)  generator  of  magnitude  4k  (newtons 

B 

or  volts  per  root  hertz).  Calculate  the  system  output  resulting  from  each  of 
these  noise  generators  separately  and  then  sum  the  squares  of  these  individual 
outputs  to  get  the  square  of  the  total  noise  output. 

Quality  Factor 

Since  mechanical  damping  is  a  critical  part  of  thermal  noise  analysis,  it 
is  useful  to  have  several  ways  of  determining  it.  In  some  cases,  the  mechani¬ 
cal  resistance  can  be  modeled  fairly  accurately  (for  example,  viscous  drag  on 
simple  shapes  at  low  Reynolds  numbers),  but  in  others  it  must  be  measured. 

For  simple  harmonic  oscillators  the  quality  factor  or  Q  is  related  to  the 
damping  and  can  often  be  measured  easily.  Several  expressions  of  Q  are 
useful : 


(1)  Q  is  2n  times  the  number  of  cycles  of  oscillation  required  for  the 
energy  of  the  oscillator  to  drop  by  the  factor  e: 

E(t)  =  exp  (-Wjjt/Q)  (5) 

(2)  Q  is  the  ratio  of  the  resonance  frequency  to  the  full  width  of  the 
resonance  peak  at  the  half-power  points  (3  dB  down  points): 


Q 


f  /Af 

0  3dB 


(6) 


(3)  Q  is  the  ratio  of  either  the  mass  reactance  or  the  stiffness  reactance 
at  resonance  to  the  resistance  in  a  series  resonant  circuit: 


Q  ■  w  m/R  ■>  k/w  R 
0  0 


(7) 


9 


NAOC-91113-50 


(4)  Q  equals  2n  times  the  energy  stored  In  the  oscillator  divided  by  the 
energy  dissipated  per  cycle.  Q  is  also  equal  to  times  the  energy 
stored  in  the  oscillator  divided  by  the  power  dissipated. 

(5)  Q  is  related  to  the  damping  factor  C  or  the  loss  tangent 

Q  =  1/(20  »  l/tan«  (8) 

Critical  damping  (R  s  Zrau^)  is  C  **  1  or  Q  ■  1/2. 

(6)  When  the  Q  is  small,  methods  (1)  and  (2)  are  not  very  useful.  In  this 
case,  the  resonance  frequency  can  be  measured  by  exciting  the  system  and 
looking  for  the  frequency  that  results  in  a  90  degree  phase  shift 
between  the  input  and  output,  ‘fhe  damping  can  be  determined  from  the 
slope  of  the  phase  change  (in  radians  per  hertz)  at  the  resonance; 


(7)  Another  approach  for  small  Q  is  to  drive  the  system  with  a  square 
wave.  The  output  waveform  will  show  "ringing"  at  the  level  transitions. 
If  the  peak- to- trough  amplitude,  a,  of  the  first  half  cycle  of  the 
ringing  and  the  trough-to-peak  amplitude,  b,  of  the  second  half  cycle 
are  measured,  then  the  damping  factor  can  be  computed  as  follows  [22, 
23]: 

In  (a  /  b) 

C  “  /—I - T 

V  In^  (a  /  b)  +  n 
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Ih£  Simple  Assg.,)>g.ram&t.9£ 

The  generic  accelerometer  sensor  Is  shown  schematically  in  Fig.  3.  The 
case  is  exposed  to  the  desired  acceleration  and  the  displacement  of  the  mass 
relative  to  the  case  (z  ■■  y  -  x)  Is  the  output  of  the  accelerometer.  Also 
shown  are  the  locr.tlon  of  the  noise  force,  and  the  free-body  diagram. 


y-x 


L  "  J 

T— T 

k(y-x)  f 

'  n 

R(y-x) 

Figure  3.  Configuration  and  free-body  diagram  for  aimpla  accelerometer. 

Of  course,  the  noise  force  also  acts  on  the  sensor  case  through  the  other  end 
of  the  damper  but  this  is  not  important  here.  The  analysis  is  simpler  in  the 
frequency  domain  where  the  signal  excitation  displacement  is  YCf),  the 
response  is  Z(f),  and  the  noise  force  is  F  .  (As  a  reminder,  y  — >  Y,  y  — » 

n 

**  2 

i«Y,  and  y  Y  in  the  frequency  domain. )  To  get  the  noise  response,  set 

the  signal  Y  to  zero  and  solve  for  in  terms  of  F^;  to  get  the  signal 
response,  set  F  to  zero  and  solve  for  Z  in  terms  of  the  case  displacement  Y. 


Solving  the  simple  accelerometer  for  the  noise  response  gives 


|Z^(f) I  «  Y4k  JR  G(f )  /  k 


(11) 
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where 

G(f )  -  1  /  /  )*]“  ♦  (f/f  )Vq* 

0  0 

and  ■■  ■<  k/m.  The  factor  G(f)  lo  plotted  In  Figure  4  for 

0  0  r  • 

4,  and  100. 


Figure  4.  Normalized  roeponee,  G(f),  of  a  Mes-apring  eaelllator 
ae  a  function  of  nonaallzed  frequency,  three  valuee  of  Q. 


In  the  accelerometer  limit  (f  «  f^),  the  noise  displacement 


|Z  (f)| 


VSITtr  /  k 

B 


which  can  also  be  expressed  in  terms  of  the  Q  of  the  oscillator: 


|Z  (f)|  »  ✓4k  T/«  kQ 

'  n  '  BO 


/ 


Ak,.T/w>Q 
B  0 


The  signal  response  is 


Z^(f)|  «  (f/f^)^  G(f)  jY^I 


(12) 

Q's  of  O.S, 


becomes 

(13) 

(14) 

(15) 


12 


IIADC-9J 113-50 


but  is  the  magnitude  of  the  Input  acceleration,  a^,  so 

|Z  (f)|  -  a  G(f)  /  (16) 

and  G(f)  Is  one  In  the  accelerometer  limit. 

The  slgnal-to-nolse  ratio  at  any  frequency  is  then 


Z  /  Z  ^  «  a\Q  /  4k  T« 

■  n  ■  BO 


(17) 


(Note:  a  moving-coil  accelerometer  has  its  electrical  output  proportional  to 
dz/dt  or  uZ;  this  is  true  for  both  signal  and  noise  however,  so  the  slgnal-to- 
noise  ratio  is  unchanged. )  Here,  the  "signal"  is  taken  to  be  the  power 
spectral  density  of  the  ambient  noise  of  the  environment.  This  is  appropriate 
for  a  sensor  designed  for  optimum  reception  of  real  signals  in  this  ambient; 
if,  however,  the  real  signals  are  guaranteed  never  to  approach  the  environmen¬ 
tal  ambient,  then  the  minimum  real  signal  level  should  be  used.  (That  could 
be  a  single-frequency  level  rather  than  a  spectral  density  in  which  case  the 
noise  power  should  be  adjusted  by  the  analysis  bandwidth. ) 

From  this  expression  for  signal-to-nolse  ratio  it  is  apparent  that  the 
SNR  can  be  improved  by: 

(1)  increasing  the  Q 

(2)  reducing 

(3)  Increasing  m  (which  also  lowers 

(4)  reducing  the  stiffness  (which  lowers  u^) 

Normally,  dropping  far  enough  so  that  it  is  within  the  band  of 
expected  signals  is  not  wise  because  this  Introduces  a  nonlinear  phase  into 
the  system  response.  Increasing  the  Q  can  cause  problems  also;  if  the 
oscillator  has  a  high  Q,  then  out-of-band  oscillations  can  be  large  (at  the 
resonance  frequency,  Q  times  the  motion  well  below  the  resonance  for  the  same 
excitation)  and  the  mechanical  system  must  have  enough  dynamic  range  to  handle 
these  large  movements.  As  will  be  discussed  later,  feedback  can  be  applied  to 
reduce  the  mechanical  range  required  but  then  feedback  circuit  must  have  an 
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adequate  dynamic  range.  Moreover,  the  feedback  circuit  Introduces  additional 
noise. 


Example  1.  Determine  the  thermal  noise  limit  for  a  critically 
damped  accelerometer  that  must  be  sensitive  enough  to  detect  the  SSO-lOdB 
levels  from  0.1  to  1000  Hz.  Assume  that  an  SNR  >  1  is  adequate.  (Unless 
otherwise  stated,  all  accelerometer'-type  hydrophones  will  be  assumed  to 
be  neutrally  buoyant  so  that  the  fluid  acceleration  equals  the  sensor 
acceleration. ) 

Ihe  lowest  practical  resonance  frequency  is  1000  Hz  and  the  Q  is  1/2 
for  critical  damping.  The  smallest  acceleration  to  be  measured  (see 
Table  I)  is  5.3  x  lO”®  m/s'*  per  i/5z,  which  occurs  at  10  Hz.  In  order  to 
have  an  SNR  >  1  at  10  Hz,  the  mass  must  be  greater  than  72  grams  (from 

Eq.  17).  Clearly,  this  would  be  a  problem  for  a  very  small  sensor. 

Critical  damping  gives  good  response  linearity  at  the  expense  of 
thermal  noise.  If  the  Q  were  raised  to  100,  the  minimum  mass  drops  to 
0.36  gm.  If  the  Q  were  kept  at  1/2  but  the  resonance  were  dropped  to  10 
Hz,  the  minimum  mass  drops  to  0.72  gm.  In  this  latter  case,  however,  the 

signal  response  now  rolls  off  above  10  Hz.  This  is  not  a  problem  for 

mechanical  thermal  noise  since  the  SNR  is  not  a  function  of  frequency, 
but  this  would  be  a  dangerous  strategy  since  other  noise  sources  could 
eventually  dominate. 

Example  2.  What  is  the  minimum  Q  required  for  a  simple  accelero¬ 
meter  with  m  =  40  X  10~*  kg  and  f^  =  300  Hz  so  that  thermal  noise  is  at 

or  below  SSO-lOdB  from  0.1  Hz  to  f  ? 

0 

Smallest  SSO-lOdB  acceleration  is  at  10  Hz  so  SNR  must  be  one  or 
greater  at  10  Hz.  From  Eq.  17,  the  Q  must  be  at  least  270.  If  the  mass 
were  reduced  to  2  x  10~^  kg,  then  the  minimum  Q  would  be  5400!  This 
illustrates  one  of  the  problems  associated  with  miniature  accelerometers. 
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Ihs  simple  Pressure  Sensor 

In  the  simple  pressure  sensor,  the  moving  mass  Is  directly  exposed  to  the 
incident  acoustic  pressure  and  the  displacement  (or  speed)  of  this  mass  is 
measured.  Such  a  sensor  is  shown  in  Fig.  5. 


Figure  S,  Schematic  diagram  of  a  elnple  preeaure  eeneor. 


For  an  arbitrary  force  F,  the  displacement  response  Z  is 


|Z]  «  F  G(f)  /  k 


(18) 


For  signal,  the  force  F  =  p^S,  where  S  is  the  area  of  the  transducer  face  and 

p^  is  the  spectral  density  of  the  signal  pressure;  for  noise,  the  force  F  ■ 

✓4k  TR,  so  the  signal-to-nolse  ratio  is 
0 


Z  /  Z  ^  «  (p  S)^  /  4k  TR 

an  a  B 


Q(p  S)^  /  4k  Tw  m 

a  BO 


(19) 


In  this  case,  the  slgnal-to-noise  ratio  can  be  improved  by: 

(1)  increasing  the  area,  A 

(2)  increasing  Q 

(3)  reducing  the  resonance  frequency 

(4)  reducing  the  mass  (unless  the  stiffness  is  also 

reduced  this  increases  u  since  w  -  k/m) 

0  0 

Reduction  of  the  mass  would  seem  to  favor  microminiaturization  but  the  u  m 

0 

product  goes  as  ✓in,  which  makes  this  dependence  weak,  whereas  the  dependence 
on  area  is  as  area  squared. 


15 


NADC'91 113-50 


Example  3.  Design  a  critically  damped  pressure  sensor  with  thermal 
noise  at  or  below  the  SSO-lOdB  curve  from  0.1  to  1000  Hz  if  the  sensor 
area  is  1  cm  . 

Here  Q  is  1/2.  Set  f  to  1000  Hz  (lowest  value  without  dropping 
resonance  into  desired  band).  The  worst  case  for  pressure  on  the 
SSO-lOdB  curve  is  at  1000  Hz  where  p^  *  53  pPa/vliz.  From  Eq.  19,  the 
mass  must  be  less  than  140  grams.  (Be  careful  with  units:  use  pressure  in 
pascals,  not  micropascals,  in  Eq.  19. )  If  the  sensor  area  were  reduced 
to  1  mm^,  then  the  mass  would  have  to  be  less  than  14  micrograms.  Also, 
the  stiffness  would  have  to  be  reduced  to  keep  the  resonance  at  1000  Hz. 

One  of  the  advantages  of  working  with  the  Q  of  a  sensor  is  that,  if  the  Q 
is  measured,  it  automatically  Includes  all  of  the  relevant  damping  terms. 

While  this  is  only  strictly  true  for  slngle-degree-of-freedom  systems  with 
frequency- independent  resistance,  many  sensors  can  be  approximately  modeled 
this  way.  In  order  for  a  Q  measurement  to  be  useful,  though,  the  Q  must  be 
measured  with  the  sensor  in  the  envlroiunent  in  which  it  will  ultimately  be 
used.  For  example,  the  Q  of  a  hydrophone  should  be  measured  in  water,  not  in 
air.  Sometimes  the  difference  may  be  negligible,  but,  especially  for  nonstan¬ 
dard  designs,  significant  loss  mechanisms  (and,  therefore,  fluctuation 
producers)  could  be  overlooked  if  the  complete  mechanical  configuration  is  not 
measured  in  situ. 

Dissipation  Mechanisms 

One  of  the  keys  to  evaluating  mechanical-thermal  noise  is  in  understand¬ 
ing  the  sources  of  dissipation  in  the  system.  In  terms  of  fluctuations,  any 
mechanism  that  allows  energy  to  escape  from  the  orderly  motion  of  the  sensor 
counts  as  dissipation.  These  mechanisms  include  mechanical  damping  in  the 
spring  and  supports,  viscous  drag,  acoustic  reradiation,  electrical  leakage, 
and  magnetic  hysteresis. 
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Since  it  Is  often  very  difficult  to  separate  various  damping  mechanisms 
in  a  device,  measurement  of  the  systom's  Q  is  a  valuable  technique  for  noise 
assessment.  As  long  as  the  damping  Is  independent  (or  nearly  so)  of  frequen¬ 
cy,  the  Q  gives  the  mechanical  resl Uance  directly  in  terms  of  the  resonance 

* 

frequency  and  either  the  mass  or  stii'fness  (Q  ■  w^m/R  «  k/w^R).  Also,  the  Q 
is  the  reciprocal  of  the  loss  tangent  or  the  reciprocal  of  twice  the  damping 
factor.  Many  times,  measurements  of  a  vibrating  system’s  dominant  mass  and  Q 
are  simple;  the  resistance  and,  therefore,  the  fluctuating  force  can  be 
calculated  from  these  measurements. 

For  capacitive  sensors,  the  very  thin  gaps  that  permit  low  polarization 
voltage  and  high  capacitance  per  unit  area  also  lead  to  squeeze-film  damping: 
the  viscous  loss  associated  vith  squeezing  the  fluid  out  between  moving 
surfaces  [24-26].  Squeeze-film  damping  can  easily  dominate  the  dissipation 
mechanisms  for  gaps  of  several  microns.  For  two  parallel  disks  of  area,  S, 
with  average  spacing,  h^,  the  equivalent  mechanical  resistance  is 

“  3pS^/2nh^  (20) 

where  p  is  the  fluid’s  viscosity  (20  x  10“®  kg/m  s  for  air  at  20*C;  lO’^  kg/m 
s  for  water  at  20*’c).  Notice  the  strong  dependence  on  the  spacing.  If  one 
disk  is  perforated,  the  damping  can  be  redviced  considerably  Il27]. 

If  a  moving  object  is  not  near  another  surface  land  the  Reynold’s  number 
is  very  low  as  would  be  expected  for  acoustic  motion),  the  mechanical  resist¬ 
ance  is  given  by  Stoke’ s  formula  (28,  29].  For  a  disk  of  radius,  a,  moving 
broadside  (a  reasonable  model  for  an  accelerometer  mass:  or  a  pressure  sensor’s 
diaphragm).  Stoke’ s  formula  gives 

R^  =  16  n  a  (21) 

where  tj  is  the  fluid’s  viscosity  (20  x  lO”®  kg/m  s  for  air  at  20*C;  kg/m 

s  for  water  at  20*0. 

Radiation  resistance  can  be  a  significant  dissipation  mechanism  above  !i 
or  10  kHz  and,  for  micronachined  sensors  with  very  small  diaphragm  trasses. 
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radiation  mass  loading  can  greatly  exceed  mechanical  mass.  The  radiation 
impedance  for  a  simple  pressure  sensor  is  approximated  closely  by  [11] 


p  c  n  a*  (ka)*  /  4  +  1  0.6  ka 


(22) 


where  k  is  the  acoustic  wave  number  (w/c),  p  and  c  are  the  density  and  sotind 
speed  of  the  fluid,  and  a  is  the  radius  of  the  piston.  (The  convention 
is  assumed.)  Notice  that  the  equivalent  resistance  (the  real  part  of  Z^)  is  a 
function  of  frequency;  consequently  the  resulting  noise  is  not  white  but, 
instead.  Increases  with  frequency.  On  the  other  hand,  the  equivalent  radia¬ 
tion  mass  loading  (the  imaginary  part  of  divided  by  w)  is  not  a  function  of 
frequency. 


The  radiation  Impedance  for  an  accelerometer  Immersed  in  a  fluid  is 
approximately  that  of  an  oscillating  sphere  [30} 

Z  »  p  c  n  a^  [(ka)*  /  3  ♦  i  2  ka  /  3]  .  (23) 

MECH  I  I 


For  any  sensor,  the  reradiation  affects  the  external  surface  of  the 
transducer  and  the  expression  for  this  equivalent  noise  pressure  does  not 
contain  any  of  the  properties  of  the  transducer.  Consequently,  this  equiva¬ 
lent  noise  can  and  usually  is  included  in  the  ambient  sea  (or  air)  noise 
pressure  [31,  32]: 


P 


rad 


rad 


/  A 


V4k_TR  '  /  A 


rad 


2  f  Vk„Tlip/c 
B 

1.8  X  10"*  f  [pPa/Vl?i] 


(24) 


in  water.  This  result  also  applies  to  immersed  accelerometers  if  they  are 
neutrally  buoyant.  The  effective  noise  pressure  is  well  below  the  SSO-lOdB 
curve  over  the  entire  range  of  0.1  to  1000  Hz.  This  nol.se  component  dominates 
ocean  ambient  noise  above  50  or  100  kHz  so  it  is  only  a  concern  for  very  high 
frequency  underwater  systems. 
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(Note:  The  same  result  for  p  can  be  obtained  by  calculating  the 
pressure  that  results  from  giving  each  normal  mode  in  the  ocean  an  energy  of 
k^T/2.  The  reradiatlon  component  is  sometimes  mistakenly  taken  to  be  the 
entire  thermal  noise  component  in  a  sensor  but  this  ignores  all  of  the 
microscopic  energy-storage  modes  in  the  sensor  structure  and  fill  fluid.  From 
the  standpoint  of  thermal  excitation,  a  microscopic  (molecular)  mode  is  no 
less  important  than  a  macroscopic  mode  like  an  acoustic  mode  in  the  ocean. ) 

A  interesting  variety  of  thermal  noise  is  that  associated  with  leaky 
dielectrics  [33].  This  dissipation  Is  characterized  by  a  loss  tangent  (tan  5) 
that  can  be  nearly  independent  of  frequency  over  a  wide  range  of  frequencies. 
The  current-voltage  relationship  for  a  capacitor  with  a  leaky  dielectric  is 

I  =  1  6>  C  V  exp(-15)  (25) 

The  real  part  of  the  equivalent  impedance  (for  small  6)  is  tanS/wC  so  the 
open-circuit  noise  voltage  is 

=  4  k  T  tan3  /  w  C  (26) 

oc  B 

The  thermal  noise  associated  with  this  dissipation  has  a  spectrum  that  is  not 
white  but,  rather,  1/f.  (Clearly,  the  loss  tangent  cannot  be  independent  of 
frequency  all  the  way  down  to  DC  because  this  would  require  infinite  fluctua¬ 
tion  power. )  Many  capacitive  micromachined  sensors  use  air  (or  vacuum)  as  the 
dielectric  where  the  loss  tangent  is  negligible;  however,  in  those  cases  where 
a  liquid  might  be  used  (to  resist  hydrostatic  pressure  in  a  hydrophone,  for 
example),  this  noise  component  should  be  calculated  especially  since  the 
micromachined  sensor  is  likely  to  have  very  small  capacitance  (V^  is 

OC 

proportional  to  1/C). 

Example  4.  A  fluid-filled  capacitive  hydrophone  has  a  receiving 
sensitivity  of  -185  dB  (IV/pPa)  and  a  capacitance  of  400  pF.  If  the  loss 
tangent  of  the  dielectric  fluid  is  0.01,  what  is  the  equivalent  noise 
pressure  corresponding  to  the  open-circuit  voltage  noise  due  to  the  leaky 
dielectric? 
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From  Eq.  26,  the  open-circuit  voltage-squared  is  6.4  x  10"^*  /  f 
volts-squared  per  hertz.  At  1  Hz,  this  is  -132  dB(V).  From  the  receiv¬ 
ing  sensitivity,  the  corresponding  pressure  is  60  dB(lpPa^/Hz) ,  which  is 
well  below  the  SSO-lOdB  curve.  This  noise  component  increases  with 
decreasing  frequency  but  the  ocean  ambient  also  increases  (faster)  with 
decreasing  frequency. 

The  following  examples  illustrate  why  Ignoring  thermal  noise  in  conven¬ 
tionally-sized  underwater  sensors  has  not  created  a  problem: 

Example  5.  A  standard  piezoceramic  hydrophone  (pressure  sensor)  has 
the  following  characteristics:  f^  =  1600  Hz,  Q  *  5,  m  «  0.4  x  10”^  kg, 
and  8^3  cm^.  If  it  is  to  be  used  from  0.1  to  1000  Hz,  will  thermal 
noise  be  a  problem? 

Again,  the  worst  case  is  at  1000  Hz.  From  Eq.  19,  the  slgnal-to- 
noise  ratio  at  1000  Hz  is  20  000  so  thermal  noise  is  not  a  problem. 

Example  6.  A  conventional  geophone  can  be  modeled  as  a  simple 
accelerometer.  If  a  geophone  123]  has  f^  =  28  Hz,  m  «  2.2  gm,  and  a 

damping  factor  <;  -  0.18,  will  thermal  noise  be  a  problem? 

*8  2 

Here,  the  worst  case  is  at  10  Hz  where  a  >  5.3  x  10  m/s  .  (Use 

S 

Eq.  8  to  calculate  Q. )  From  Eq.  17,  the  slgnal-to-noise  ratio  at  10  Hz 
is  6.  While  this  is  not  as  large  as  in  the  previous  example,  thermal 
noise  would  not  be  a  problem  in  this  case  either. 

Mechanical  thermal  noise  is  often  overlooked  as  limiting  noise  source 
because  conventionally  sized  underwater  sensors  have  thermal  noise  well  below 
the  quietest  background  noise  in  the  ocean.  Thermal  noise  can  be  a  serious 
issue  though  for  very  small  sensors  as  Examples  1  through  3  illustrate.  Small 
physical  size  should  not  be  the  criterion  for  deciding  whether  or  not  to 
consider  mechanical-thermal  noise;  this  noise  component  should  be  evaluated 
for  any  sensor  that  is  intended  for  extreme  sensitivity.  For  example, 
mechanical-thermal  noise  is  routinely  considered  in  the  design  of  gravity-wave 
detectors  [34]  even  though  these  detectors  are  physically  very  large. 
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Noise  Estimation  from  Sensor  Frequency  Response 

Sometimes  the  information  available  for  a  sensor  Is  insufficient  to 
directly  calculate  Q  or  R.  One  consequence  of  equipartition  is  that  the 
details  of  damping  need  not  be  knovm.  The  total  thermal  energy  is  fixed  by 
equipartition;  the  system's  frequency  response  can  only  redistribute  that 
energy.  If  a  frequency  response  curve  (measured  or  theoretical)  is  given  for 
the  sensor,  this  curve  can  be  used  to  evaluate  the  mechanical-thermal  noise 
for  simple  sensors. 

In  order  for  this  to  work,  the  spectral  shape  of  the  signal  response  of 
the  sensor  must  be  the  same  as  the  spectral  shape  of  the  output  of  the  sensor 
due  to  thermal  agitation  of  the  sense  element.  For  the  simple  accelerometer, 
the  spectral  shape  of  the  acceleration  response  is  identical  to  the  shape  of 
the  noise  response  (Eq.  11);  for  the  simple  pressure  sensor,  the  shape  of  the 
pressure  response  is  identical  to  the  shape  of  the  noise  reponse  (Eq.  18). 

This  may  not  be  true  for  more  complicated  sensors  but,  if  the  behavior  is 
dominated  by  a  single  mass-spring  system,  then  the  following  procedure  may  be 
suitable  as  a  first  approximation  to  the  noise  characteristics. 

The  spectral  density  of  the  noise  is  generally  not  constant  with  frequen¬ 
cy  because  the  mechanical  spring-mass  system  shapes  the  spectrum  but  the 
broadband  noise  is  still  governed  by  equipartition  so 

00 

<  x^  >  -  S  x^(f)  df  -  k  T  /  k  (27) 

B 

0 

(Note:  this  is  the  appropriate  starting  point  if  the  sensor  output  is  proport¬ 
ional  to  the  displacement  of  the  sensor  mass;  if,  as  in  the  case  of  a 
moving-coil  accelerometer,  the  output  is  proportional  to  the  velocity  of  the 
mass,  then  <v  >  and  Eq.  3  should  be  used  to  write  the  velocity  equivalent  of 
Eq.  27. )  This  expression  says  that  the  energy  in  the  mode  associated  with  the 
mean-square  displacment  is  constant  regardless  of  its  spectral  distribution. 
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Figure  6  shows  the  spectral  distribution  uf  the  mean-square  displacement 
for  three  systems  with  the  same  k  (so  the  right-hand  side  of  Eq.  27  remains 
constant),  but  different  resonance  frequency  or  Q.  The  spectral  shape  is 
different  in  each  case  but  the  area  under  each  curve  is  the  same.  For 
example,  increasing  the  Q  does  not  reduce  the  total  noise  energy,  but  it 
concentrates  the  energy  near  the  resonance  thereby  reducing  it  in  the  band 
below  resonance. 
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Figure  6.  Nolee  power  dtetrlbutlon  Tor  three  oectlletorc  with  the  eenr 
spring  constant  but  different  resonance  frequencies  and  Q’s. 


If  the  given  sensor  power  response  function  is  called  G(f)  and  the 
displacement  noise  has  the  same  shape,  then  x^(f)  ■  c^G(f)  and  all  that  is 
needed  is  the  value  of  the  constant  c^.  This  is  found  directly  from  Eq.  27: 
the  integral  is  simply  the  area,  A^,  under  the  response  curve  so 

c.  -  k  T  /  k  A,  (28) 

1  B  f 


and 
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x^(f) 

■  G(f )  k  T  /  k  A 

B  r 

(When  the  area  Is  calculated,  the  curve  must  be  expressed  In  a  power-llke 
quantity  —  displacement  squared,  pressure  squared  —  and  certainly  never  in 
dB. )  Once  the  displacement  noise  is  known,  the  factor  relating  sense  mass 
displacement  to  sensor  output  can  be  used  to  calculate  the  sensor  output 
noise.  For  a  moving-coil  sensor  whose  output  is  proportional  to  the  velocity 
of  the  moving  element,  the  equivalent  velocity  noise  expression  is 

v^(f)  =  G(f)  k  T  /  m  A  (30) 

B  r 

Example  7.  A  pressure  sensor  with  an  equivalent  spring  constant  of 
10^  Nt/m  has  a  flat  frequency  response  between  0  and  500  Hz  and  rolls  off 
rapidly  beyond  that.  Evaluate  its  mechanical-thermal  noise. 

The  response  can  be  approximated  by  G(f)  1  from  0  to  500  Hz  and 
G(f)  •=  0  above  500  Hz.  The  area  under  the  curve  is  then  500  and  Eq.  29 
then  gives  x  ■  9  x  lO”^®  m/VHz.  If  the  active  face  of  the  sennor  were  3 
mm  by  3  mm,  then  the  equivalent  noise  pressure  would  be  the  equivalent 
noise  force,  kx,  divided  by  this  area.  In  this  case,  the  noise  pressure 
is  56  fiPa/vliz,  which  is  roughly  equal  to  the  SSO-lOdB  curve  at  500  Hz. 


Sensors  with  Feedback 

One  of  the  ways  of  reducing  the  mechanical-thermal  noise  is  to  increase 
the  Q  of  the  sensor.  As  has  been  discussed,  a  high  Q  can  cause  dynamic  range 
problems  in  that  out-of-band  signals  can  produce  large  excursions  of  the  sense 
mass  and  these  excursions  must  be  accomodated  by  the  mechanical  design  in 
order  that  the  in-band  signals  not  be  distorted. 
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There  is  a  means  for  eliminating  this  mechanical  problem  and  that  is  by 
negative  feedback:  the  output  of  the  sensor  is  amplified  and  applied  to  u 
force-generating  mechanism  that  opposes  the  motion  of  the  sense  element.  If 
properly  done,  the  feedback  force  keeps  the  sense  element  virtually  stationary 
and  the  sensor  output  is  taken  from  the  input  to  the  feedback  force  generator. 
This  is  sometimes  knovm  as  force-balance  feedback. 

In  this  way,  the  mechanical  frequency  response  is  made  effectively  flat 
and  the  mechanical  dynamic  range  problem  is  eliminated.  Even  if  the  sensor 
output  without  feedback  is  not  linearly  related  to  the  sense  mass  motion  (the 
electron-tunneling  sensor  is  inherently  nonlinear),  the  sensor  output  with 
feedback  is  linear. 

These  improvements  are  not,  however,  free.  If  the  phase  response  of  the 
feedback  loop  is  not  carefully  controlled,  the  intended  negative  feedback  can 
become  positive  at  some  frequencies  and  the  sensor  can  be  driven  into  self 
oscillation.  Also,  negative  feedback  does  not  reduce  the  effects  of  thermal 
noise;  in  fact,  the  feedback  loop  adds  additional  noise.  Finally,  while  the 
mechanical  dynamic  range  problem  is  eliminated,  the  feedback  circuitry  must 
have  sufficient  electrical  dynamic  range  to  successfully  force-balance  the 
mass  motion  through  the  resonance.  Fortunately,  electrical  dynamic  range  is 
usually  easier  to  come  by  than  mechanical  dynamic  range. 

The  implementation  of  feedback  in  the  simple  accelerometer  is  shown 
schematically  in  Fig.  7. 


Figure  7.  Simple  accelerometer  with  feedback  force  applied. 
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The  force,  f ,  is  applied  to  oppose  the  motion  of  the  mass.  (Two  ways  In  which 
the  force  can  be  applied  are  electrostatically  or  by  means  of  a  plezoceramlc 
element.)  The  control  circuit  (using  frequency  domain  notation)  would  look 
something  like  the  circuit  In  Fig.  8. 


Figure  8.  Control  circuit  for  f orce-balence  feedback. 

A  is  the  gain  of  a  very  high  gain  dlsplacement-to-voltage  conversion  stage  and 
B  is  the  transfer  constant  for  the  voltage-to**force  transducer. 

In  the  frequency  domain,  the  response  Z  to  a  case  displacement  X  Is 
Zg  =  w^m  X  /  [^  i/m  -  IwR  -  (k  +  AB)  J  (31) 

which,  for  very  large  A  (i.e. ,  AB  »  «^m),  becomes 

Z  «  w^m  X  /  AB  (32) 

s 

Since  A  Is  very  large,  Z^  is  kept  very  small  by  the  feedback  force.  The 
actual  sensor  ou'.put  is  the  voltage  input  to  the  feedback  force  transducer. 
This  voltage  Is 

V  «  X/B  »  ma/B  (33) 

out  • 

in  terms  of  thf.t  case  acceleration,  a  .  Notice  that  the  response  to  signal  is 
flat  as  long  as  the  frequency  Is  low  enough  so  that  AB  >>  (i>^m;  the  higher  the 
feedback  gain,  the  higher  the  high-frequency  rolloff  is.  Besides  removing  the 
response  peak  around  the  resonance,  the  high  frequency  rolloff  can  be  delayed 
by  high  gain  in  the  feedback  loop. 
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For  noise  responfie,  the  sense  mass  displacement  Z  is 

/  AB  (34) 

and  the  corresponding  voltage  output  is 

V  ,  «  V^^JR  /  B  (35) 

out  B 

Consequently,  the  i;lgnal~to-noise  ratio  Is  obtained  by  squaring  Eqs.  33  and  35 
and  dividing; 


V. 


2 


2  2 
i.  m 


4k  TR 
s 


2 

a^mQ 

4kju 

B  0 


(36) 


which  is  identical  to  the  expression  without  feedback  giver  in  Eq.  17.  So 
feedback  dees  not  improve  the  sensor’s  immunity  to  then  a L  noise.  In  fact, 
components  in  the  feedback  circuit  will  introduce  additUnal  noise.  Because 
feedback  permits  use  of  high  Q,  the  introduction  of  addil  ional  noise  in  the 
feedback  loop  may  be  more  than  offset  by  the  noise  reduction  from  increased  Q. 


Compound  Sensors 


The  foregoing  analyses  apply  to  sensors  that  are  well-represented  by  a 
single  mass-spring  system.  While  this  is  adequate  for  many  hydrophones,  some 
designs  are  more  complex.  The  solution  procedure  for  more  complicated  systems 
is  more  tedious  but  it  is,  in  principle,  identical  to  the  simpler  systems. 

For  example,  consider  an  accelerometer  with  two  mass-spring  components. 

A  schematic  diagram  and  the  relevant  fr<ie-l,jdy  diagrams  are  showit  in  Fig.  9. 


j 
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Figure  9.  Configuration  end  freo-bcdy  dlagrens  for  a 
compound  accelerometer. 


As  before,  the  analysis  is  more  conveniently  done  In  the  frequency  domain.  To 
get  the  signal  response,  set  F  .  and  F  to  zero  and  solve  for  S  (  *  X  ••  Z) 

nl  nZ  V 

as  a  function  of  the  case  displacement  X.  To  get  the  noise  response,  set  X  to 


zero  and  set  F  to  zero  and  solve  for  S  as  a  function  of  F  ;  then  set  X 
nZ  nl  nl 


and  F  to  zero  and  solve  for  S  as  a  function  of  F  .  Since  noise  powers 

nl  ^  n^  _  n2 

add,  the  total  noise  is  =  S®  +  .  Also,  F  »  VAk  TR  and  F 

n  nl  n2  nl  B  1  n2 

✓4k  TR  . 

B  2 


In  the  accelerometer  limit  (w  «  the  signal  response  is 
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and  the  slgnal-to-nolse  ratio  is 


SNR 


a^m  0  (1  +  m./m.  +«Vw*)^ 

■  11 _ 2  1  12 _ 

4k  Tu  1  +  (m  /in,)(Q/Q,)(w/M,)® 

Bl^  121212 


(38) 


If  m  <<  m  .  u  »  u  ,  and  the  k‘s  and  the  Q*s  are  about  equal  for  the  two 
2  12  1 

mass-spring  systems,  the  factors  in  square  brackets  are  approximately  equal  to 
one  and  the  expressions  reduce  to  those  of  a  single  mass-spring  system. 


Complex  mechanical  devices  are,  perhaps,  best  analyzed  by  first  drawing 
the  electrical  equivalent  circuit  and  then  using  a  computer-based  circuit 
analysis  program  to  determine  the  frequency  and  noise  response.  This  approach 
is  discussed  In  Appendix  B. 


Limits  of  Anpllcabllltv 


Just  because  a  sensor  passes  the  thermal  noise  test  does  not  mean  that  It 
will  not  have  a  noise  problem.  One  of  the  other  noise  factors  discussed  In 
subsequent  sections  may  dominate.  However,  If  It  fails  the  thermal  noise 
evaluation,  then  the  sensor  does  have  a  problem. 

The  noise  terms  are  not  really  uniformly  distributed  over  all  frequencies 
—  that  would  Imply  Infinite  power.  More  correctly,  the  factor  k  T  should  be 

B 

replaced  by  [13,  17] 


_ hf _ 

exp  (hf/k.T)  -  1 

where  h  Is  Plank’s  constant  (6.6  x  10”^*  J  s).  For  hf  «  k^T  this  expres¬ 
sion  reduces  to  k  T  but  above  hf  ■  k  T  In  frequency,  the  noise  power  Is 
reduced.  At  room  temperature,  this  rolloff  occurs  above  about  10*^  Hz;  so, 
this  is  is  not  a  concern  for  typical  underwater  acoustic  sensors. 
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Some  authors  [35-37]  Include  the  zero-point  energy  term  in  Eq.  39,  which 
amounts  to  an  additional  hf/2.  Whether  or  not  this  term  should  be  included  is 

13 

the  subject  of  some  debate  [37]  and,  since  it  is  not  significant  below  10 
Hz,  it  is  only  of  academic  interest  here.  While  its  Inclusion  does  produce 
the  "proper"  limit  as  T  goes  to  zero  it  also  makes  Eq.  39  diverge  at  high 
frequency.  The  zero-point  energy  does  put  an  ultimate  limit  on  optical 
systems  and  cryogenlcally  cooled  systems  [38,  39]  but  the  zero-point  energy 
cannot  be  extracted  from  the  system  [35]  and  so,  from  the  standpoint  of  noise 
power  delivered  to  the  observer,  need  be  considered  only  as  a  fundamental 
limit  on  the  observability  of  a  process. 


Shot  Noise 


If  a  system  is  in  equilibrium,  then  thermal  noise  provides  a  complete 
description  of  its  Internal  fluctuations.  There  are,  however,  several 
important  non-equilibrium  sources  of  noise  [18];  One  of  these  is  shot  noise. 


If  a  signal  is  carried  by  discrete  "particles"  (electrons  or  photons,  for 
example)  and  these  particles  are  emitted  randomly,  then  there  will  be  noise 
associated  with  that  signal.  In  the  case  of  an  electrical  current,  the 
spectral  density  of  the  noise  component  of  that  current  is  [40,  41] 


1  =  ✓  2  I  e  [amps/Vliz] 

N 


(40) 


where  I  is  the  average  current  flowing  and  e  is  the  charge  on  an  electron  (1.6 
X  10”^®  coulombs  per  electron).  This  is  shot  noise. 


Photons  generated  by  a  laser,  charge  carriers  crossing  potential  barriers 
at  semiconductor  Junctions,  and  electrons  emitted  from  the  filament  in  a 
vacuum  tube  are  all  examples  of  random  emissions  that  lead  to  shot  noise.  On 
the  other  hand,  the  current  produced  by  applying  a  voltage  to  a  conductor  is 
spatially  correlated  over  large  distances  and  results  in  very  little  shot 
noise  [19,  40,  42,  43].  Shot  noise  results  from  charge  carriers  crossing 
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potential  barriers  of  some  sort  and  the  process  of  Jumping  the  barrier  Is 
random. 

In  a  bipolar  Junction  transistor,  the  charge  carriers  are  randomly 
injected  into  the  emitter-base  and  base-collector  depletion  layers  and  it  Is 
in  these  layers  that  the  carriers  interact  with  externally  applied  voltages 
[18],  These  layers  are  too  thin  (i.e. ,  the  transit  times  of  the  carriers  are 
too  short)  to  permit  the  carriers  to  reach  equilibrium  with  the  "stationary" 
molecules.  Consequently,  the  noise  produced  in  the  external  circuits  by  these 
Junctions  is  directly  related  to  the  randomness  of  the  transit  of  carriers 
across  the  Junction  rather  than  being  related  to  thermal  vibrations  in  the 
semiconductor  lattice.  This  noise  is  primarily  shot  noise. 

In  a  field-effect  transistor  (FET),  the  carriers  interact  with  the 
external  voltages  in  a  channel  that  is  long  enough  to  allow  the  carriers  to 
reach  thermal  equilibrium  with  the  channel.  Consequently,  the  characteristics 
associated  with  the  random  emission  into  the  channel  are  lost  and  the  noise  in 
a  field-effect  transistor  is  primarily  thermal  noise. 

Even  if  the  photon  generation  by  a  laser  were  not  random,  detection  of 
that  light  by  a  photodiode  involves  random  emission  of  electrons  in  response 
to  the  incident  photons,  so  the  photodiode  output  has  shot  noise.  A  typical 
value  for  photodiode  sensitivity  is  one  microamp  output  for  two  microwatts  of 
light  input:  this  value  is  not  arbitrary  but  is  the  result  of  each  incoming 
photon  (with  a  frequency  in  the  visible-light  range)  forcing  a  single  electron 
out. 


Example  8.  Suppose  a  fiber-optic  interferometric  accelerometer  [44] 
uses  200  fiU  of  optical  power,  which  produces  100  pA  of  current  out  of  the 
photodiode.  The  mass  is  0.S4  kg,  the  resonance  frequency  is  240  Hz,  and 
the  Q  is  10.  Interferometer  sensitivities  are  often  given  in  terms  of 
radians  of  phase  shift  (in  the  interference  fringes)  per  sensed  quantity. 
In  this  case,  the  phase  shift  per  unit  displacement  of  the  sense  mass  is 
8500  radians  per  micrometer  of  displacement.  The  sensitivity  can  be 
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increased  simply  by  adding  more  fiber  to  arms  of  the  interferometer.  If 
the  laser  noise  is  50  pradians  per  Vtiz,  is  there  any  point  to  increasing 
the  sensitivity? 

This  sensor  has  both  shot  and  thermal  noise  components.  Analyze  the 
shot  noise  first.  An  interference  fringe  goes  from  light  to  dark  so 
there  would  be  (at  most)  a  lOr  |iA  current  change  corresponding  to  n 
radians  (u  half  cycle}  of  fri.ige  shift.  'The  electrical  sensitivity  is 
then  lOO/TT  or  about  30  /iA/rad.  The  worst-case  shot  noise  is  for  an 
average  current  of  100  pA  so  the  shot  noise  spectral  density  (Eq.  40)  is 
i  *  6  X  10~^^  A/Vfiz  and  the  minimum  detectable  signal,  in  terms  of 

n 

fringe  shift,  is  then  .  »  1  /(30  pA/rad)  »  0.2  prad/^/Hz.  This  is  well 

below  the  laser  noise  so  this  sensor  is  not  shot-noise  limited. 

The  noise  displacement  of  the  sense  mass  resulting  from  thermal 
noise  can  be  calculated  from  Eq.  13.  The  displacement  noise  is  9. 3  x 
lO"^®  m/VHz.  For  an  optical  sensitivity  of  8500  rad/pm,  the  correspond¬ 
ing  phase  noise  is  B  prad/'v^.  Although  this  is  greater  than  the  shot 
noise,  the  laser  nolfie  st'.ll  dominates.  Consequently,  the  overa}.l 
sensitivity  of  this  accelerometer  coulu  be  increased  (with  more  fiber) 
somewhat. 

Example  9.  Accelerometers  are  not  sensitive  to  hydrostatic  pressure 
(that  is,  until  the  case  .Implodes!)  but  a  pressure  sensor  is  unless  it  is 
somehow  compensated  for  the  hydrostatic  pressure.  Another  way  of  saying 
this  is  that  the  simple  pressure  sensor  responds  all  the  way  down  to  DC 
pressure.  Suppose  that  an  uncompensated  fiber-optic  pressure  sensor  has 
been  designed  to  operate  between  100  and  300  meters  depth  over  the 
frequency  range  0.1  to  100  Hz.  Suppose  that  this  is  not  an  interfero- 
mef:rlc  sensor  but  that  the  variation  in  light  output  can  be  adjusted  so 
that  it  ranges  smoothly  from  dark  to  bright  over  the  range  of  pressure 
sensed.  Evaluate  the  shot  noise  performance  [45] . 

Since  the  sensor  is  not  compensated  for  hydrostatic  pressure 
Vc.rlations,  it  must  respond  to  pressures  over  the  entire  hydrostatic 
renge  —  about  2  MPa  hei'e.  (Uach  UO  meters  in  depth  adds  about  1  MPa  to 
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the  hydrostatic  pressure.)  If  is  the  maximum  current  out  o'  the 
detecting  photodiode,  then  the  sensitivity  Isj 

M  -  I  /  AP  (41) 

pd 

where  AP  Is  the  total  required  range  of' pressure  to  be  sensed  (2  MPa  In 

this  case).  I  Is  also  the  current  from  which  to  compute  the  (worst- 
pd 

case)  shot  noise  spectral  density  1^  from  Eq.  40,  so  the  minimum 
detectable  pressiure  Is 

p  s  1  /  M  =  AP  '  2  e  /  I  (42) 

nln  n  pd 

The  worst  case  In  the  band  of  Interest  Is  at  100  Hz  where  the  required 
p  is  200  pPa/VHz  so  the  required  photodiode  current,  I^,  would  be  32 
amps,  which  would  require  an  64  watt  laser!  In  order  for  this  sensor  to 
be  a  practical  underwater  sensor,  some  means  would  have  to  be  developed 
to  eliminate  the  DC  pressure  response. 

The  shot  noise  was  particularly  serious  In  this  last  example  because  the 
sensor  responded  to  pressure  all  the  way  to  zero  frequency  (hydrostatic  or  DC 
pressure).  A  conventional  air-backed  piezoceramic  pressure  hydrophone  does 
not  respond  to  DC  pressure  because  leakage  currents  In  the  sensor  remove  the 
charge  produced  by  very  slow  straining  of  the  material.  Also,  the  plezoceram- 
Ic  sensor  does  not  modulate  shot-noise  producing  current,  so  shot  noise  Is 
not  a  consideration. 

The  spectrum  of  shot  noise  Is  uniformly  distributed  in  frequency 
("white")  up  to  rather  high  frequencies.  In  order  for  the  shot  noise  relation 
to  hold,  there  must  be  many  charge  carriers  per  cycle.  There  are  I/e  elect¬ 
rons  per  second  in  a  current,  I,  or. there  are  I/ef  electrons  per  cycle  at  the 
frequency,  f.  For  a  current  of  1  nA,  there  is  one  electron  per  cycle  at  6  x 
lo’  Hz,  so  for  frequencies  well  below  this,  the  shot  noise  spectrum  would  be 
white.  Increasing  the  current  increases  the  high-frequency  limit. 
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1/f  Noise 

Except  for  a  very  high  frequency  cutoff,  the  fundamental  mechanisms  of 
both  thermal  and  shot  noise  are  uniformly  distributed  in  frequency.  There  is. 
however,  a  very  commonly  occuring  noise  component  whose  spectral  distribution 
is  not  flat  but,  instead,  drops  with  increasing  frequency  so  that  the  power 
distribution  is  approximately  f~*.  This  sort  of  behavior  is  seen  in  measure¬ 
ments  of  the  base  current  of  transistors,  the  speed  of  ocean  currents,  the 
flow  of  sand  in  an  hourglass,  the  yearly  flow  of  the  Nile  over  the  iast  2000 
years,  traffic  flow  on  expressways,  sunspot  activity,  and  the  loudn;;ss  of 
classical  music  [40,  46,  47]. 

Unfortunately,  there  is  no  simple  physical  model  capable  of  explaining 
this  noise  component  as  there  is  in  the  case  of  thermal  or  shot  noJse  [48, 

49] .  The  fact  that  this  1/f  distribution  of  power  occurs  so  often  in  so  many 
vastly  different  settings  suggests  that  the  underlyls\g  explanation  should  not 
depend  on  the  specific  details  of  any  single  physical  process.  In  addition, 
the  1/f  portion  of  the  spectral  distribution  can  extend  over  many  lecades  so 
the  time  scales  required  to  explain  this  behavior  must  range  from  milliseconds 
to  hours  or  even  days  in  some  instances  [47]. 

The  distribution  is  not  precisely  f”^  either.  Observations  of  different 
processes  show  that  the  exponent  can  range  [49]  from  -0.8  to  -1.4.  Further¬ 
more,  when  the  exponent  is  -1  or  more  negative,  some  low-frequency  cutoff  must 
be  present  or  there  would  be  infinite  power  in  the  spectrum  [33,  48,  49]. 

One  attempt  to  model  the  spectral  shape  of  this  noise  component  relies  on 
a  distribution  of  relaxation  times  for  the  process  [49].  Roughly  speaking, 
relaxation  time  is  the  time  it  takes  a  system  to  return  to  “normal"  after 
being  disturbed.  If  the  system  reverts  to  its  undisturbed  state  immediately 
(i.e. ,  it  has  no  “memory"),  then  the  associated  power  spectrum  is  independent 
of  frequency.  If  the  system  has  one  characteristic  relaxation  time,  x,  then 
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the  spectral  shape  is 

w(f)  «  A  T  /  II  +  (wt)®]  (43) 

where  A  is  an  arbitrary  constant.  This  spectrum  varies  as  1/f^  for  ut  »  1 
(or  f  »  l/Znx) . 

If,  for  some  reason,  a  system  has  a  distribution  of  relaxation  times  and 
the  probability  distribution  of  those  t's  has  Just  the  right  shape,  then  a  1/f 
spectrum  results.  For  example,  if  the  relaxation  time  is  controlled  by  some 
sort  of  energy,  E,  t  is  likely  to  be  exponentially  related  to  that  energy  (t  » 
T  exp(E/k  T));  if  the  probability  distribution  of  the  energy  is  uniform  over 

0  B 

some  range  of  t,  then  the  spectrum  will  have  1/f  slope  over  the  equivalent 
frequency  range.  A  reasonable  physical  model  for  1/f  noise  in  semiconductors 
has  been  developed  based  on  such  a  distribution  of  relaxation  times  where  the 
relaxation  processes  are  connected  with  trapping  and  release  of  the  charge 
carriers  ("generation-recombination"  noise)  [41,  49].  This  approach  has  not 
been  so  successful  in  describing,  for  example,  1/f  noise  in  metallic 
conductors. 

In  electrical  devices,  the  1/f  noise  is  proportional  to  the  current 
passing  through  the  device.  Some  authors  have  suggested  that  the  noise  results 
from  random  fluctuations  in  resistance  that  are  sensed  by  the  current  rather 
than  from  fluctuations  caused  by  the  current  [33,  48,  49].  This  conclusion 
contradicts  the  notion  that  1/f  noise  is  a  non-equilibrium  phenomenon;  some 
external  source  of  energy  would  be  required  to  keep  the  resistor  from  equilib¬ 
rium  and  the  most  likely  suspect  is  the  current.  In  some  materials,  the  noise 
is  proportional  to  the  volume  of  material,  but  the  noise  is  often  sensitive  to 
surface  conditions  as  well. 

Any  component  that  carries  current  is  a  potential  source  of  1/f  noise 
[18]  in  addition  to  thermal  noise  [37].  Carbon  resistors  produce  much  more 
1/f  noise  than  wirewound  resistors  and  so  should  be  avoided  (in  favor  of 
either  wirewound  or  metal-film  resistors)  in  low-noise  preamplifiers.  Contact 
resistance  —  bad  solder  Joints,  terminal  lugs  screwed  to  aluminum  frames  — 
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can  contribute  significant  amounts  of  1/f  noise.  Poor  selection  of  capacitor 
types  for  coupling  or  bypass  capacitors  (especially  through  leakage  currents) 
can  contribute  also:  electrolytic  capacitors  are  particularly  poor;  ceramic 
capacitors  are  moderately  poor;  poly  ‘ster  are  acceptable;  and  teflon,  poly¬ 
styrene,  and  polypropylene  capacitors  are  excellent.  Polarized  capacitors 
should  be  avoided  where  possible.  For  applications  requiring  large  capaci¬ 
tance,  tantalum  capacitors  are  acceptable.  Any  polarized  capacitor  must  be 
protected  against  reverse  DC  bias,  which  may  occur  during  power-on  transients, 
because  they  can  emit  noise  pulses  up  to  several  hours  afterward  [50]. 

Because  there  is  no  simple  predictive  technique  to  estimate  the  magnitude 
of  1/f  noise,  its  effects  must  generally  be  measured.  Fortunately,  manufact¬ 
urers  often  provide  sufficient  information  about  transistors  and  amplifier 
chips  to  allow  estimation  of  noise  (of  all  sources  including  1/f).  The  next 
section  summarizes  how  to  use  this  information.  Before  turning  to  that  topic 
though,  it  is  worth  pointing  out  that  1/f  noise  is  not  only  a  low  frequency 
problem.  Oscillators  must  be  nonlinear  in  order  for  their  amplitude  to  be 
stable  and  this  nonlinearity  coupled  with  1/f  noise  results  in  sidebands  on 
the  oscillator  output  [18]. 

It  is  also  possible  that  1/f  noise  is  not  limited  to  the  electrical  part 
of  a  sensor.  Mechanical  parts  in  the  sensor  assembly  may  also  contribute  1/f 
noise  according  to  one  study  of  condensor  microphones  [51-53].  In  this  case, 
the  energy* required  to  keep  the  system  from  equilibrium  was  likely  to  have 
been  the  polarizing  voltage. 

An  effective  way  of  avoiding  the  effects  of  1/f  noise  in  very-low- 
frequency  measurement  systems  is  to  arrange  the  signal  to  modulate  a  high- 
frequency  reference  signal.  For  example,  a  capacitive  sensor  can  be  used  in 
an  AC  bridge  circuit.  If  the  sensor  is  in  one  leg  of  the  bridge  and  a 
reference  capacitor  is  in  the  other,  then  the  output  is  a  differential  signal: 
oscillator  noise  effectively  cancels  out  since  its  signal  flows  through  both 
branches. 
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Even  better  If  the  sensor  itself  can  be  designed  with  a  "push-pull" 
output.  A  capacitive  accelerometer  could  be  built  with  the  moving  plate 
suspended  between  two  stationary  plates.  The  two  capacitors  so  formed  would 
be  used  In  the  bridge  circuit,  one  In  each  leg. 

Another  approach  Is  to  use  an  Inductor  to  form  a  tuned  circuit  that  is 
driven  electrically  at  Its  resonance.  The  phase  shift  resulting  from  capaci¬ 
tance  changes  Is  directly  proportional  to  the  electrical  Q  (see  Eq.  9).  This 
approach  does  not  by  itself  remove  the  oscillator  noise  but  a  bridge  circuit 
can  be  constructed  In  which  each  leg  of  the  bridge  contains  a  tuned  circuit 
[43]. 


Amplifier  Noise 

Once  the  sensing  element  itself  has  been  designed  for  adequate  noise 
performance,  It  is  usually  necessary  to  amplify  the  resulting  signal  to  a 
useable  level.  While  electrical  noise  calculations  are  much  more  familiar 
than  mechanical  noise  calculations,  It  Is  worth  summarizing  the  principles 
here  because  a  poor  preamplifier  can  negate  the  noise  performance  of  a  quiet 
sensor. 

The  safest  way  to  analyze  the  noise  performance  of  a  preamplifier  circuit 
connected  to  a  sensor  Is  to  calculate  the  slgnal-to-nolse  ratio  directly  from 
the  manufacturer’s  data  (or  measurements)  for  the  equivalent  Input  noise 
voltage  and  current  for  the  transistors  or  chips  [18,  40,  54].  This  Is  the 
procedure  that  will  be  outlined  below,  but  some  more  common  but  less  useful 
definitions  will  be  mentioned  first  In  order  to  warn  of  some  dangers  and  to 
make  connections  between  alternate  expressions  of  noise  performance. 

One  of  the  most  commonly  quoted  noise  measures  Is  noise  figure.  Noise 
figure  Is  ten  times  the  common  logarithm  of  the  noise  factor;  noise  factor  for 
a  particular  device  Is  the  ratio  of  the  signal-to-nolse  power  at  the  Input  of 
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the  device  to  the  signal-to-noise  power  at  the  output  of  the  device: 


noise  factor  »  nf  -  SNR  /  SNR 

INPUT  OUTPUT 

noise  figure  ■  NF  ■  10  log  (nf) 


(44) 


Defined  this  way,  noise  factor  Is  always  greater  than  one  (and  noise  figure  Is 
always  positive)  and  Is  a  measure  of  the  degradation  In  slgnal-to-nolse  ratio 
caused  by  the  device. 


This  sounds  like  a  good  approach  but,  In  practice,  device  manufacturers 
often  quote  only  the  optimum  noise  figure.  Noise  figure  Is  a  strong  function 
of  source  resistance,  so  If  the  source  resistance  Is  not  given  for  that  noise 
figure  or  the  stated  source  resistance  Is  not  close  to  the  actual  source 
resistance  used,  the  noise  figure  Is  of  little  value. 


Equivalent  InpMi  112iS£  g9Virg?s 


Of  far  more  value  are  the  manufacturer’s  specifications  for  equivalent 
Input  voltage  and  current  noise  for  the  device.  By  convention,  these  values 
are  equivalent  values  at  the  input  to  the  device  and  so  are  added  to  the 
source  noise  contributions  (thermal,  shot,  1/f).  The  output  levels  can  be 
found  simply  by  multiplying  by  the  device  gain  but  the  slgnal-to-noise  ratio 
remains  the  same.  In  effect,  the  real  device  has  been  replaced  by  a  noise- 
free  amplifier  and  voltage  and  current  noise  generators  as  In  Fig.  10. 


REAL  LEV  ICE 


NOISE-FREE 

DEVICE 


Figure  10.  Input  noise  current  and  voltage  generators. 
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Because  the  noise  voltage  and  current  are  given  as  the  equivalent  input 
values,  the  noise  figure  can  be  determined  from  the  input  noise  alone  (the 
ratio  of  the  input  signal  to  the  input  signal  is  one!).  If  the  device  were 
truly  noise-free  (i  and  e  equal  to  zero),  then  the  input  noise  voltage 

n  n 

squared  would  be  4k  TR  from  the  source  resistance,  R  .  In  terms  of  the  noise 

B  ■  '  ■ 

figure,  the  actual  input  noise  for  the  real  device  would  then  be 

e^  *  (4k  TR  )  nf  (45) 

B  ■ 

or  (dB//V*).  ,  *=  10  log  (4k  TO  )  +  NF 

Input  not  so  Bn 

This  can  be  misleading,  though,  because  the  input  noise  may  not  be 
entirely  due  to  the  obvious  input  circuit  resistance.  The  input  noise  also 
includes  the  mechanical  thermal  noise  of  the  sensor  and,  maybe,  shot  noise  and 
1/f  noise.  Of  course,  an  equivalent  source  resistance  could  be  derived  to 
include  all  these  effects  but  this  is  an  awkward  approach.  Again,  the  best 
approach  is  to  calculate  the  equivalent  input  noise  directly  and  either  forget 
about  noise  figure  altogether  or  calculate  the  noise  figure  directly  from  the 
equivalent  input  noise. 

If  the  source  noise  can  be  represented  by  an  equivalent  resistance,  R^, 
then  the  amplifier  can  be  analyzed  using  the  model  of  Fig.  11. 


R 


Figure  11.  Hodel  for  nolso  analyts  of  amplifier. 

(Note:  if  the  source  impedance  is  complex,  take  the  real  part  for  R  for  those 

s 

terms  represent Ing  generation  of  thermal  noise  —  4k  TO  terms  —  but  use  the 

B 

magnitude  of  the  source  impedance  for  those  terms  representing  voltages 
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generated  by  current  flowing  through  the  source  Impedance  —  I  R  terms,  for 

n 

example.)  The  equivalent  input  noise  voltage  and  current  (40,  5S~S7]  for 
several  low-noise  op  amps  and  one  low-noise  preamp  (LM381)  are  shown  in 
Figures  12  and  13. 


Frequency  (Hz) 


riQure  12.  Input  vultaqa  nal*«  for  Intcgrnt.ad  circuit*, 

line*  Indicate  dope  correepondlng  to  l/f  nolae. 
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For  preamplifier  chips  and  transistors  without  feedback,  the  equivalent 
input  noise  voltage  resulting  from  the  amplifier  is 

e^  =  +  (i  R  )*  (46) 

n  n  « 

(Sometimes  it  is  more  convenient  to  calculate  the  amplifier’s  input  noise 
current 


i^  «  (e  /R  )^  +  )  (47) 

n  a  n 

If  R  represents  the  entire  noise  contribution  of  the  source,  then  the  total 
input  noise  voltage  is  e^  from  Eq.  46  plus  4k^TR^.  This  total  noise  voltage 
es  a  function  of  source  resistance  is  shown  in  Figure  14  for  the  devices  of 
Figures  12  and  13. 
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Notice  that,  in  each  case,  there  is  an  optimum  source  resistance  that 
results  in  the  smallest  noise  excess  over  source  resistor  noise.  Notice  also 
that  this  optimum  resistance  is  substantially  different  from  device  to  device. 
For  example,  the  LT1028  is  particularly  good  for  source  resistances  between 
100  and  1000  ohms  while  the  LF411  is  good  for  source  resistances  between  1  and 
10  megohms.  For  best  noise  performance,  it  is  not  enough  to  select  a  low' 
noise  chip;  it  is  also  necessary  to  match  the  chip  to  the  source  Impedance. 

Operat^on^;  ampliO^rfi  WUh  LggdbagJs 

The  curves  in  Figure  14  represent  the  noise  performance  for  the  op  amps 
without  feedback.  This  is  not  realistic  since  op  amps  are  almost  always  used 
with  feedback  and  the  feedback  circuits  add  additional  noise  (both  from 
thermal  noise  in  the  feedback  resistors  and  interaction  of  the  amplifier  i^ 
and  e  with  the  feedback  resistors).  For  the  noninverting  feedback 

n 
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configuration  (see  Appendix  C)  shovm  in  Fig.  15, 


Figure  IS.  Nonlnverttng  configuration  for  on  operational  anpllflor. 

the  voltage  gain  Is  G  *  1  total  Input  nolso  voltage  (including 

source  resistor  noise)  Is 

e^  «*  +  4k  T(R  +  R  )  +  l^CR  +  R  (48) 

n  Bop  nop 

where  R  Is  the  parallel  combination  of  R  and  R  :  R  ■  R  RV(R  +R  ). 

p  12pl212 

For  the  Inverting  configuration  of  Fig.  16, 


Figure  16.  Operational  anpltfier  In  Inverting  configuration. 


the  voltage  gain  G  «  ^2^1  total  input  noise  voltage  is 

e*  -  (1  +  R  /R,)*  e®  +  4k  TO 

a2n  Ba  (49) 

+  4k  TR,(R  /R  )*  ♦  iV 

B  2  a  2  n  a 

where  R  is  the  sum  of  R  and  R  .  A  more  complete  discussion  of  noise 

a  si 

analysis  of  op  amps  with  feedback  and  JFET  amplifiers  is  in  Appendix  C. 
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By  comparing  Eqs.  48  and  49  to  Eq.  46.  it  is  obvious  that  feedback 
increases  the  amplifier  noise  contribution.  Once  the  amplifier  configuration 
has  been  chosen,  either  Eq.  47  or  48  can  be  used  to  calculate  the  input  noise 
to  which  must  be  added  any  sensor  noise  component  not  included  in  the  source 

resistance  R  .  If  the  sensor  noise  component  is  an  equivalent  current  noise, 

•  2 

the  amplifier  noise  may  be  more  conveniently  expressed  as  a  current:  i  * 

e^/R^.  In  any  case,  the  proper  R^  may  also  include  mechanical  resistance. 

The  best  way  to  insure  that  all  components  are  accounted  for  is  to  draw  the 
complete  electro-mechanical  equivalent  circuit  as  described  in  Appendix  D. 

(Note:  Some  low-noise  preamplifier  chips,  such  as  the  LM381  shown  in 
Figures  12-14,  have  optional  differential  Inputs.  If  the  differential  Inputs 
are  used,  the  e  value  for  the  chip  must  be  doubled  (i  stays  the  same).) 

n  n 

Example  10.  A  non-traditional  way  of  sensing  the  displacement  of 
the  sense  mass  in  an  accelerometer  is  by  electron  tunneling.  If  a  probe 
is  placed  very  close  (10  ^  to  10  meters)  to  a  conducting  surface  and  a 

small  voltage  is  applied  from  the  conducting  surface  to  the  probe,  a 
current  flows  because  of  wave-like  tunneling  of  electrons  across  the  gap 
[9,  58].  The  actual  mechanism  of  current  generation  is  unimportant  here; 
the  Important  detail  is  that  a  current  is  generated  that  is  proportional 
to  displacement.  (The  fact  that  the  proportionality  is  not  linear  will 
be  ignored. ) 

Consider  an  electron-tunneling  accelerometer  with  m  =  25  x  10"*  kg, 
a  resonance  frequency  of  200  Hz  and  a  Q  of  200.  The  change  in  tunneling 
current  with  a  change  in  displacement  is  7  amps  per  meter  and  the  average 
tunneling  current  is  1  nanoamp  with  an  applied  voltage  across  the  gap  of 

7 

100  mV.  The  preamplifier  is  an  inverting  op  amp  with  R^  «  0  and  R^  ■  10 
ohms.  (In  this  configuration,  the  op  amp  acts  as  a  current- to-voltage 
convertor  with  an  output  voltage  equal  to  10  times  the  input  current. ) 
Evaluate  the  noise  performance  and  select  a  suitable  op  amp. 

Since  this  is  a  current-producing  sensor,  it  is  convenient  to 
express  all  noise  components  in  equivalent  currents.  The  mechanical- 
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thermal  noise  displacement  is  given  by  Eq.  13,  and,  when  multiplied  by 
the  displacement-to-current  factor  of  7  A/m,  the  resulting  noise  current 
is  2.8  X  lO"^^  A/vliz.  The  shot-nolse  current  from  Eq.  39  is  l.S  x  10“** 
A/^IE.  The  tunneling  gap  has  an  equivalent  resistance  equal  to  the  bias 
voltage  (100  mV)  divided  by  the  tunneling  current  (1  nA)  or  10  ohms; 
however,  this  contributes  no  noise  because  tunneling  is  not  an  equilibrium 
process:  there  is  no  mechanism  that  can  bring  the  tunneling  waves  into 
equilibrium  as  they  transit  the  gap.  Without  considering  the  preamplifier 
contribution  then,  this  device  is  limited  by  mechanical-thermal  noise. 

(The  mechanical-thermal  contribution  is  slightly  larger  than  the  equivalent 
SSO-lOdB  value  at  10  Hz  but  that  problem  will  be  ignored  in  this  example.) 

The  preamplifier  noise  can  be  determined  from  Eq.  49  divided  by  the 
source  resistance  squared  (to  get  i  )  and  with  »  0: 


i 


2 


1^  +  4k  T/R 

n  B  2 

+  (1  +  R  /R)^/  R^ 

n  ■  2  m 


(50) 


The  first  and  third  terms  on  the  right-hand  side  reflect  the  impact  of  the 
amplifier  current  and  voltage  noise  components.  The  second  term  results 
from  thermal  noise  in  the  feedback  resistor  and  the  third  term  includes 
the  effects  of  the  amplifier  voltage  noise  on  the  feedback  resistor. 

In  order  to  illustrate  the  selection  problem  with  cp  amps,  three 
types  will  be  considered:  two  expensive  chips  —  LT1028  and  AD549,  and 
one  cheap  chip  —  the  LF411.  The  equivalent  noise  voltage  and  current 
for  each  chip  can  be  estimated  from  Figures  12  and  13.  At  10  Hz,  the 
values  are  given  in  Table  II. 

Table  II.  10  Hz  no  lee  values  for  LT102B,  AD549,  and  Lr411. 


LT1028 

AD549 

LF411 

e  Iv/^Hz] 

1  X  10"’ 

1  X  lO""^ 

6  X  10'® 

i  [A/vTIz] 

5  X  10"^^ 

1  X  10”*® 

2  X  10‘^* 

n 
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Table  III  sununarizss  the  Individual  current  noise  components  (in 
amps^Hz  so  that  they  Ciin  be  added  directly)  for  the  entire  system  with 
each  of  the  three  choices  of  op  amp.  The  terms  on  the  right-hand  side  of 
Eq.  50  have  been  named  amp.Ufier  current,  feedback  thermal,  and  amplifier 
voltage  In  order  of  their  appearance  in  that  equation. 


Table  III.  Current  noliis  coaponenta  for  aanaor/aepl It'!  ar  syatea. 


LT1028  AD549  LF411 

7.8  :<  10‘^*  7.8  x  10‘^®  7.8  x  10"“ 

3.2  X  lO"^®  3.2  X  10"^®  3.2  X  10“®® 

1.6  X  10"^'  1.6  X  lO"*"^  1.6  X  10"^ 

1.2  X  10“®^  1.2  X  10"^®  4.4  X  10"®’ 

2.5  X  10“®®  1.0  X  10"®®  4.0  X  10"®® 

For  each  op  amp,  the  dominant  noise  component  is  boxed.  Clearly,  the 
expensive  LT1028  would  be  an  extremely  poor  choice  since  its  noise 
component  is  much  higher  than  the  intrinsic  sensor  noise.  While  the 
LT1028  does  have  exceptionally  low  voltage  noise,  its  current  iioise  is 
quite  high.  Since  the  source  resistance  is  very  high  in  this  applica¬ 
tion,  it  is  necessary  to  select  a  chip  with  low  current  noise.  (The 
LT1028  would  be  excellent  for  applications  with  very  low  source  resist¬ 
ance.  )  Both  the  ADS49  and  the  LF411  would  perform  well  since  the  sensor 
noise  dominates  in  both  cases  but  the  LF411  is  considerably  cheaper. 

This  analysis  has  been  done  for  10  Hz,  which  is  the  worst-case 
frequency  for  designing  to  the  SSO-lOdB  curve  for  the  simpL^  accelero¬ 
meter  alone.  The  preamplifier  introduces  1/f  noise,  however,  so  the 
lowest  frequency  of  interest  should  also  be  examined.  If,  for  example, 
the  system  Is  intended  to  be  used  down  o  0.1  Hz,  the  amplifier  voltage 
and  current  noise  would  be  conslderabl'<  higher.  If  the  manufrcturer  does 
not  supply  equivalent  noise  va!tues  at  a  low  enough  frequency,  the  values 
can  be  estimated  by  extrapolating  the  given  values  by  1/f  (in  power). 

The  voltage-squared  (or  current-squared)  noise  would  be  100  times  greater 
at  0.1  Hz  than  at  10  Hz.  This  would  increase  the  values  in  the  last  two 
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lines  of  the  above  table  by  a  factor  of  100.  The  LF411  and  the  ADS49 
systems  would  still  be  adequate. 

General  Considerations  £&£  Pr^MPUfUrg  [18.  40.  SO.  54] 

(1)  The  first  stage  of  amplification  usually  dominates  the  noise  perfor¬ 
mance  of  the  electronics.  It  should  be  a  simple  preamplifier  that  passes  the 
frequency  range  of  Interest.  Mixers,  complex  filters  and  detectors  should  be 
saved  for  later  stages. 

(2)  Don't  select  an  amplifier  by  noise  figure  alone.  Dc  the  complete 
analysis  as  in  the  previous  example  and  account  for  the  actual  source 
resistance. 


(3)  Don't  add  a  resistor  In  series  with  a  source  to  more  closely  match 
the  optimum  source  resistance  for  a  particular  amplifier.  This  may  make  the 
noise  figure  of  the  amplifier  alone  lower  but  it  will  add  more  than  enough 
additional  thermal  noise  to  compensate. 

(4)  With  op  amps,  don't  forget  to  Include  the  feedback  components  In  the 
noise  calculations.  (Use  Eq.  48  or  49,  not  Eq.  46. ) 


(5)  If  the  manufacturer’s  specifications  on  e  and  1  don’t  go  low  enough 

2  2 

in  frequency,  estimate  the  appropriate  values  by  extrapolating  e  and  1  as 

n  n 

l/f.  (Do  this  even  If  the  manufacturer’s  curve  shows  a  constant  value  wli  h  no 
low-frequency  increase;  the  l/f  breakpoint  in  the  curve  Is  probably  Just  holow 
the  given  portion.)  Of  course,  don't  try  to  extrapolate  beyond  the  pa8sbai'.d 
of  the  device. 


(6)  For  systems  that  require  very  low  noise  electronics,  better  perfor¬ 
mance  can  be  obtained  with  a  discrete  (FET  or  bipolar  transistor)  input  stage 
added  tu  an  op  amp.  Here,  a  JFET  would  be  selected  for  a  hlgh-lmpedance 
source,  while  a  bipolar  transistor  would  be  selected  fcr  a  low-impedance 
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source.  MOSFCTs  are  not  used  for  low  frequency  applications  because  they  have 
large  1/f  noise  components. 


:  .  '^5;  *  '  - 


Sunnary 

As  the  proceeding  examples  illustrate,  mechanical-thermal  noise  can  have 
observable  consequences  in  micromachined  sensors.  The  gross  effects  can  be 
estimated  easily  either  through  Nyquist’s  Relation  or  the  Equipartltion 
Theorem.  Frequently,  this  estimate  is  sufficient  to  determine  if  there  is 
thermal  noise  problem.  Of  course.  Just  because  a  sensor  passes  the  thermal 
noise  test  does  not  mean  that  it  will  not  have  a  noise  problem.  There  are 
many  other  sources  of  noise  in  an  electromechanical  transducer  system; 
however,  if  the  sensor  fails  the  thermal  noise  evaluation,  then  it  does  have  a 
problem. 

If  the  sensor  can  be  approximated  as  a  simple  pressure  sensor  or  a  simple 
accelerometer,  then  an  estimate  for  the  sensor  noise  can  be  made  from  the 
sensor  Q,  from  a  calculation  of  the  dominant  damping  mechanism,  or  from  a 
measured  or  calculated  frequency  response.  Many  times  this  simple  estimate 
will  be  sufficient.  If  not,  then  a  more  detailed  mechanical  analysis  or  a 
complete  electromechanical  equivalent  circuit  may  be  required.  The  effects  of 
the  sensor  preamplifier  can  be  included  if  the  complete  circuit  is  drawn. 

In  short,  when  a  high-sensitivity  sensor  is  being  designed,  an  analysis 
of  mechanical-thermal  noise  should  be  included  at  an  early  stage  to  avoid 
being  trapped  with  an  unacceptably  high  noise  floor. 
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Appendix  A.  Derivation  of  Nyquist’s  Relation 

Nyquist’s  relation  gives  the  connection  between  the  spectral  density  of 
the  fluctuations  of  a  system  in  thermal  equilibrium  and  the  dissipation  of 
that  system.  Since  this  relation  comes  from  equilibrium  thermodynamics,  the 
result  is  independent  of  the  physical  model  used.  The  following  derivation  is 
done  for  a  damped  mass-spring  oscillator  [59-61]  but  it  can  be  done  Just  as 
well  for  a  damped  inductor-capacitor  oscillator  [17,  19]. 


As  discussed  for  the  simple  pressure  sensor,  the  displacement  response, 

Z,  for  a  damped  harmonic  oscillator  to  which  an  arbitrary  force,  F,  Is  applied 
is  given  by  Eq.  13  rewritten  here 

|Z|  =  F  /  1:  y  (1-fi^)^  +  nW  (A-1) 

where  (1  =»  w/w  . 

0 

From  equlpartition,  Eq.  2,  k<2^>  *  k  T.  The  quantity  <z^>  is  equal  to 
the  integral  of  |Z|  over  all  frequencies  so  the  equlpartition  expression  can 
be  written 


k 


df _ 

(1-n^)^  +  nW 


(A-2) 


where  F  is  the  noise  driving  force,  which  is  assumed  to  be  uniformly  distri- 
n 

buted  over  all  frequencies.  Changing  the  integration  variable  from  f  to  (2 
yields 


F*« 
n  0 

2nk 


dR 


(1-n^)^  + 


V 


(A-3) 


The  denominator  of  the  Integral  is  quadratic  in  and  so  can  be  factored. 
Then  the  integral  can  be  solved  by  complex  contour  integration  but  it  also 
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appears  in  Gradshteyn  and  Ryzhlk  [62,  Eq.  3.264.2]  where  its  value  is  fo\md  to 
be  nQ/2.  Consequently, 


n  Q 

n  0  _ 

2nk  2 


4  k 

0 


(A-4} 


or 


F  “  i^4k  TR  [force/^/Sz] 

n  B 

which  applies  even  if  R  is  a  function  of  frequency  [16]. 


(A-5) 


As  mentioned  at  the  beginning  of  this  Appendix,  the  sane  result  can  be 
obtained  by  analyzing  a  variety  of  different  physical  models.  Alternative 
systems  leading  to  the  Nyquist  relation,  Eqn.  A-S,  are 

(1)  A  spring-damper  system  acting  on  an  object  with  negligible  mass  [36, 
63].  (For  example,  Langevin’s  method.) 

(2)  A  transmission  line  with  matched  resistive  terminations  [21,  43,  64]. 
This  is  the  device  used  by  Nyquist  to  obtain  the  electrical  form  of 
his  relation  (for  "Johnson"  noise). 

(3)  Fermi-Dirac  conduction  electrons  in  a  metal  [17,  37]. 

(4)  Equilibrium  between  a  small  surface  and  an  infinite  volume  of  fluid 
[32,  65].  Tills  is  one  of  the  most  interesting  examples  because  the 
problem  can  be  solved  in  two  ways.  The  spectral  density  of  the 
pressure  fluctuations  can  be  derived  by  the  acoustic  analog  of 
blackbody  radiation  (as  is  done  for  the  Rayleigh- Jeans  Law)  where  the 
surface  is  assumed  to  be  inside  a  volume,  the  frequency  distribution 
for  the  normal  inodes  in  the  volume  is  obtained,  the  nodes  are  each 
given  their  thermal  equipartlon  energy,  and  the  volume  is  allowed  to 
increase  indefinitely.  The  result  is  [31] 

dp^  »  (p/’nk^T/c)  f*  df  (A-6) 
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The  Nyquist  relation 

dF^  =  dp^  »  4k  TR  df  (A-7) 

B 

where  S  is  the  surface  element  area,  Implies  that 
the  associated  resistance  Is 

R  =  pnS^fVc  (A-8) 

From  strictly  hydrodynamic  considerations,  the 
radiation  resistance  can  be  calculated  for  a  surface 
element  in  an  infinite  fluid  volume  and  the  result 
is  identical  to  Eq.  A~8. 
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Appendix  B.  Equivalent-Circuit  Modeling  with  SPICE 

Noiee  analysis  can  be  performed  on  complicated  mechanical  systems  easily 
by  using  a  computer  simulation  of  the  equivalent  electrical  circuit.  This 
appendix  summarizes  the  use  of  one  of  the  popular  software  tools  —  SPICE  — 
to  do  these  analyses  [66]. 

The  first  step  is  to  select  an  appropriate  analogy  between  electrical 
quantities  and  mechanical  quantities  [6S1.  For  an  entirely  mechanical  system 
the  choice  between  the  impedance  or  the  mobility  analog  is  mostly  one  of 
preference.  If,  however,  it  is  desired  to  analyze  a  mechanical  transducer 
with  its  associated  electronics,  the  natural  association  of  pressure  with 
voltage  in  a  piezoelectric  or  capacitive  sensor  favors  the  impedance  analog, 
while  the  moving-coil  sensor  favors  the  mobility  analog.  The  impedance  analog 
will  be  used  here;  relations  for  the  mobility  analog  can  be  written  similarly. 

The  Impedance  analog  is  as  follows: 

velocity,  V  < - >  current,  i 

force,  F  * - >  voltage,  e 

mass:  F  *  m  ^  < - >  ®  ®  L  ^  (B-1) 

fit  at 

(where  v  is  velocity  in  inertial  frame) 

spring:  F  =  k  Jv  dt  < - *  e  =  i  Ji  dt  (B-2) 

(where  v  is  the  velocity  difference  between  the  ends) 

damper:  F  =  R  v  i - *  e  »  R  i  (B-3) 

(where  v  is  the  velocity  difference  between  the  ends) 

From  these  relations  it  is  clear  that  Inductance,  L,  is  the  analog  of  mass,  m; 
the  inverse  of  capacitance,  1/C,  is  the  analog  of  stiffness,  k;  and  the 
mechanical  resistance,  R,  is  the  analog  of  electrical  resistance,  R. 
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The  key  to  the  Impedance  analogy  lies  in  the  distinction  betwffen  the 
relevant  velocity  for  a  mass  and  the  relevant  velocity  for  either  a  spring  or 
a  damper.  This  distinction  is  Illustrated  in  Fig.  B-1  (mechanical  symbol  on 
left,  electrical  symbol  with  appropriate  "current”  on  right). 
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Figure  B-1.  Electronechenical  Impedance  analogjr. 


When  drawing  the  electrical  circuit,  split  the  current  as  required  to  keep  the 
appropriate  velocity  or  velocity  difference  on  the  objects.  Figure  B-2 
illustrates  this  procedure  for  a  simple  example. 


Figure  B-2.  Conatruction  of  Impedance  analogy  and  creation 
of  proper  branch  currents  (velocities). 
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In  SPICE,  voltage  and  current  (corresponding  to  force  and  velocity, 
respectively)  can  be  monitored  easily.  V/hen  acceleration  and  displacement 
must  be  measured,  a  small  subterfuge  is  required.  To  measure  acceleration, 
insert  a  small  inductor,  m^,  and  measure  the  voltage  (force)  across  it.  Since 
f  >  m^a,  the  acceleration  at  that  point  is  equal  to  F/m^.  Obviously  the 
inuerted  Inauctor  must  be  small  with  respect  to  the  other  com]xinents  so  as  not 
to  influence  the  circuit  behavior. 

To  measure  displacement,  insert  a  small  capacitor,  1/k^,  and  measure  the 

voltage  (force)  across  it.  The  force  F  •  k^x  so  the  displacement  is  equal  to 

F/k  .  If  the  circuit  values  can  be  arranged  so  that  m  and  k^  can  be  set  to 
0  0  0 

one,  then  the  acceleration  and  displacement  can  be  road  directly  from  the 
voltage  across  the  probe  component. 

There  is,  however,  an  unfortunate  problem  with  SPICE:  all  nodes  must  have 
a  DC  path  to  ground.  Insertion  of  the  displacement  probing  capacitor  will 
often  Isolate  a  node  from  ground.  To  re-establish  DC  contact  with  ground,  add 
a  vary  large  inductor  in  parallel  with  the  probe  capacitor  (see  Fig.  B-3), 


I - 1 

Figure  B-3.  Use  of  Inductor  to  provide  DC  path  sround  capacitor. 

Don’t  use  a  parallel  resistor  because  that  may  upset  the  noise  calculation. 
Also,  don’t  use  a  small  test  resistor  to  measure  velocity  (current)  as  this 
can  also  upset  the  rolse  calculation.  To  measure  velocity,  insert  an  Indepen¬ 
dent.  voltage  source  with  zero  amplitude.  SPICE  allows  the  syntax  I(VTEST)  to 
display  the  current  flowing  through  the  source,  VTEST. 

The  total  voltage  noise  Is  calculated  using  the  .NOISE  statement. 
Unfortunately,  there  Is  no  convenient  way  to  calculate  the  total  current  noise 
directly. 
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Voltage  and  current  sources  in  SPICE  are  used  to  drive  the  mechanical 
system.  Pressure  Is  force  divided  by  area  (If  It  Is  uniform  on  the  transducer 
face)  so  a  voltage  source  with  an  amplitude  equal  to  the  desired  pressure 
times  area  represents  the  driving  pressure.  A  velocity  driver  Is  directly 
represented  by  a  current  source.  There  Is  no  acceleration  (or  displacement) 
source  so  a  velocity  source  should  be  used  with  a  series-connected  probe 
Inductor  (or  capacitor)  to  measure  the  acceleration  (or  displacement)  drive 
level. 

Before  considering  some  examples  of  electrical  analogs  of  mechanical 
devices,  the  acoustic  analogy  should  be  mentioned  since  some  transducers  are 
more  easily  represented  by  acoustic  lumped  parameters  rather  than  mechanical 
lumped  parameters  (111.  In  the  acoustic  analogy,  current  Is  the  analog  of 
volume  velocity  and  voltage  Is  the  analog  of  pressure.  Volume  velocity,  q,  Is 
the  total  volume  flow  rate  through  some  reference  surface  and  is  defined  for  a 
vector  velocity  field,  v, 

q  =  f  v  •  dS  (B-4) 

aria 

where  dS  is  the  differential  surface-area  element  expressed  as  a  vector  ncrmal 
to  the  surface.  Normally,  the  acoustic  analogy  is  only  used  if  the  fluid 
velocity  can  be  considered  to  be  constant  over  and  normal  to  the  area  In 
question.  In  this  case,  the  volume  velocity  is  simply  q  ®  vS  where  S  is  the 
area.  The  acoustic  analogy  is  used  for  devices  In  which  fluid  flow  is  a 
significant  property. 

Acoustic  Impedance  xn  defined  as  pressure  divided  by  volume  vfiluclty  and 
is  equal  to  the  mechanical  Impedance  (force  divided  by  velocity)  d!.vided  by 
the  area  squared.  Therefore,  mechanical  elements  are  represented  xs  follows: 

2 

mass - »  Inductor  with  L  «  m/S 

R  - »  resistor  with  R  ■  R  /S^ 

RRch  aoch 

stiffness  - 4  capacitor  with  1/C  »  k/S^ 
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In  fluid-filled  systems,  stiffness  elements  ere  often  volumes  of  fluid. 

For  gases,  these  volumes  are  represented  as  [11] 

volume  - >  capacitor  with  C  ■  V/yp^ 

where  V  is  the  volume,  y  is  the  ratio  of  specific  heats  (1.4  for  air),  end  p^ 
is  the  ambient  pressure.  For  liquids,  replace  yp^  with  the  bulk  modulus, 

B.  The  principles  for  equivalent  circuit  construction  for  acoustic  analogs 
are  similar  to  those  for  mechanical  impedance  analogs. 

Example  Simple  Pressure  Sensor.  The  mechanical-impedance 
equivalent  circuit  is  shown  in  Fig.  B-4. 


.  F  m  . .  1/k 


/////////////// 

Figure  B-4.  Equivalent  circuit  for  elnple  preeeure  eeneor. 

and  the  electrical  equivalent  circuit  with  a  displacement  probe  to 
measure  Z  (and  numbered  nodes)  is  shown  in  Fig.  B-5. 


VSRC 


Flgiure  B-S.  SMCE  equivalent  circuit  for  pressure  sensor. 


To  analyze  a  specific  case,  the  actual  mechanical  component  values 
can  be  used  and  CTEST  can  be  selected  so  that  CTEST  »  C.  In  this 
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illustration,  however,  the  values  will  all  be  normalized  as  follows: 

(1)  Set  CTEST  «•  1  so  that  the  voltage  V(4,0)  is  directly  equal  to  the 
displacement,  z. 

(2)  Pick  C  «  CTEST  (e.g.,  lE-5) 

(3)  Pick  L  so  that  the  resonant  fre^ency,  is  one.  ■  1/LC  so 

L  ■  2533.03  for  this  example.) 

(4)  Pick  R  to  set  Q:  Q  *•  u^L/R.  For  Q  ■  10  in  this  example,  R  ■ 
1591. 5!>. 

AC  analysis  gives  the  sensor  output  displacement  for  a  constant  pressure 
input;  NOISE  analysis  gives  the  thermal  noise  output  resulting  from  R. 

The  following  input  file  for  SPICE  performs  these  analyses  on  this 
circuit: 

SIMPLE  PRESSURE  SENSOR 
VSRC  1  0  AC  10 

L  1  2  2533.03 

C  23  1E~5 

R  34  1591.55 

CTEST  0  4  1 

LTEST  0  4  1E12 

.AC  DEC  10  0.01  100 

.NOISE  V(4,0)  VSRC  20 

.PRINT  AC  VM(0,4) 

.PRINT  NOISE  ONOISE(M) 

.END 

(A  SPICE  manual  should  be  consulted  for  a  detailed  explanation  of  these 
statements  [66]). 

For  a  frequency  well  below. the  resonance  (e.g. ,  f  ■  0.01),  the 
signal  output  VM(0,4)  •  lE-5  which  equals  F/k  (the  amplitude  of  the 
driving  force  is  one  here)  as  it  should.  Also,  the  noise  output 
ONOISE(M)  -  5. 137E-14  m/VRz  which  equals  V4k  TO  /  k  as  it  should.  At  the 

B 

resonance,  both  of  these  values  should  be  Q  (10)  times  larger  and  they 
are  in  the  SPICE  simulation. 
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Example  B-2.  Simple  Accelerometer.  The  mechanical -Impedance 
equivalent  circuit  Is  shovm  In  Fig.  B-6 


Figure  B-6.  Equlvelent  circuit  for  simple  eocelerometer. 


and  the  electrical  equivalent  circuit  with  an  input  acceleration  probe 
(LCASE)  and  an  output  displacement  probe  (CTEST)  Is  shown  In  Fig.  B-7. 


Figure  B-7.  SPICE  equlvelent  circuit  for  eccelerometer. 

For  normalized  output  as  In  the  previous  example,  set  CTEST  ■  1  and 
LCASE  B  1.  (Note;  Physically,  LCASE  represents  a  small  mass  associated 
with  the  accelerometer  case  and  CTEST  represents  a  very  weak  spring 
connected  between  the  mass  pointer  and  the  case  pointer.)  The  other 
component  selections  are  Identical  to  those  for  the  simple  pressure 
sensor . 

AC  analysis  gives  the  output  displacement  and  the  input  accelera¬ 
tion.  Note  that,  while  the  velocity  driver  Is  constant,  the  Input 
acceleration  Is  a  function  of  frequency  so  the  transfer  function  must  be 
calculated  by  dividing  the  output  displacement  V(0,4)  by  the  Input 
acceleration  V(2,l).  For  the  pressure  sensor,  the  Input  force  was 
constant  (and  equal  to  one)  so  the  transfer  function  (the  sensor’s 
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receiving  sensitivity)  between  Input  pressure  and  output  dlsplaceaent  was 
equal  to  the  output  displacement  times  the  transducer  area.  The  SPICE 
NOISE  analysis  gives  the  thermal-noise  displacement  resulting  from  R. 

The  following  Input  file  for  SPICE  performs  these  analyses  on  this 
circuit: 

SIMPLE  ACCELEROMETER 

ISRC  0  1  AC  1  0 

L  02  2533.03 

C  2  3  lE-5 

R  34  1S9J.SS 

LCASE  1  2  1 

CTETT  0  4  1 

LBYPS  0  4  1C12 

.AC  DEC  10  0.01  100 

.NOISE  V(4,0)  ISRC  20 

.PRINT  AC  VM(1,2)  VM(0,4) 

.PRINT  NOISE  ONOISE(M) 

.END 


(The  magnitude  of  the  displacement-output  to  acceleration-input  ratio  Is 
VM(0,4)/VM(1,2). )  The  values  for  the  sensor  signal  response  and  noise 
output  from  the  SPICE  simulation  agree  with  the  theoretical  calculations. 
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Example  0*3.  Compound  Ucceleroxmter.  The  mechenlcel’-iapedence 
equivalent  circuit  Is  shown  In  Fig.  B-8 


and  the  electrical  equivalent  circuit  Is  shown  In  B-9. 


LCASr 


I5RC 


Figure  B-9.  SPICE  equivalent  circuit  for  compound  acceleroMter. 


The  following  input  file  for  SPICE  performs  the  signal  and  noise 
analysis  on  this  compound  accelerometer.  The  values  were  selected  so 
that,  for  the  oscillator,  the  resonance  frequency  Is  one  and  the 

Q  Is  5;  and,  for  the  oscillator,  the  resonance  fre^ency  Is  two 

and  the  Q  Is  2.5. 
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COMPOUND  ACCELEROMETER 
ISRC  0  1  AC  1  0 

LI  i  3  2533.03 
Cl  2  4  lE-5 
R1  46  3183.10 

L2  30  1266.51 

C2  3  5  lE-5 
R2  56  3183.10 

LCASE  1  2  1 

CTEST  6  0  1 

LBYPS  6  0  1E12 

.AC  DEC  10  0.01  100 

.NOISE  V(6.0)  ISRC  20 
.PRINT  AC  VM(1,2)  VM(6,0) 

.PRINT  KOISE  ONOISE(M) 

.END 


Th«  displacement  (signal)  response  and  the  noise  spectral  density 
for  this  compound  system  as  calculated  by  the  SPICE  simulation  are  shown 
in  Fig.  B-10.  Notice  that  the  displacement  response  and  the  noise 
spectral  density  do  not  have  the  same  dependence  on  frequency.  This 
would  introduce  an  error  (although  small  in  this  case)  in  any  estimate  of 
noise  power  based  on  a  measured  dlsplaicement  response  curve. 
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Response 


Normalized  Frequency 


Figure  B-10.  Signal  and  nolae  reaponsc  coaputad  by  SPICE  For 
compound  accalaronator. 


Much  more  elaborate  mechanical  systems  can  be  modeled  with  SPICE.  With 
some  ingenuity  other  forms  of  friction  and  even  backlash  can  be  described  with 
various  combinations  of  SPICE  elements  (67). 


I 


(Times  10E-13) 
Noise  Spectral  Density 
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Appendix  C.  Noise  Anslysls  of  Op  Anps  and  JFETs 

There  are  several  subtle  problems  that  arise  In  calculating  the  noise  of 
combined  sensor/ampllf ier  systems.  In  this  Appendix,  the  noise  characteris¬ 
tics  of  op  amps  and  JFET  preamplifiers  will  be  considered  in  sufficient  detail 
to  present  these  problems  and  provide  a  model  for  analysis  of  similar 
electronic  circuits. 

The  usual  representation  of  noise  in  an  amplifier  or  a  transistor  is  as 
Independent  voltage  and  current  sources.  When  the  device  is  in  a  circuit, 
especially  if  feedback  is  applied,  it  may  not  be  possible  to  maintain  this 
Independent  voltage  and  current  representation.  (For  example,  the  total  noise 
current  would  contain  the  amplifier  noise  current,  but  the  total  noise  voltage 
could  also  include  a  term  resulting  from  the  amplifier  noise  current  flowing 
through  a  circuit  resistor. )  Fortunately,  it  is  not  necessary  to  do  this.  In 
fact,  it  is  better  to  analyze  the  circuit  with  the  signal  source  connected  and 
calculate  the  total  noise  voltage  (or  the  total  noise  current)  without  trying 
to  artificially  develop  separate  voltage  and  current  components. 

Two  of  the  more  useful  techniques  for  performing  this  noise  analysis  are: 

(1)  Compute  the  output  voltage  resulting  from  each  noise  source  (amplifier 
current,  amplifier  voltage,  resistor,  shot)  separately;  divide  by  the 
circuit  gain  to  refer  the  values  back  to  the  input;  and,  sum  the  squares 
of  the  individual  contributions  to  get  the  total  mean-square  input  noise 
voltage. 

(2)  Use  the  siginal-to-noise  ratio  theorem  for  feedback  circuits  [43]:  The 
output  signal -to-noise  ratio  is  independent  of  the  load  impedance.  ITie 
output  SNR  based  on  output  signal  and  noise  currents  with  the  output 
shorted  to  ground  is  identical  to  the  output  SNR  based  on  voltages  with 
the  output  open-circuited.  Shorting  oi  open-circuiting  the  output  can 
simplify  the  calculations  considerably. 

To  Illustrate  these  techniques,  the  noninverting  op  amp  circuit  will  be 
analyzed  by  the  first  procedure  and  the  inverting  op  amp  will  be  analyzed  by 
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the  second  procedure.  For  these  analyses,  the  simple  model  for  an  op  amp 
(Fig.  C-1)  will  be  used  [40j. 


Figure  C-1.  Slaple  «odel  for  operetlone!  eapllfler. 

In  this  model,  A  Is  very  large.  (When  the  output  Is  short-circuited  to  ground 
In  the  second  procedure,  a  sniall  series  resistor  representing  the  open-loop 
output  Impedance  of  the  op  anp  will  be  Inserted  to  keep  the  current  finite. ) 

Noninverting  Amplifier 

The  noninverting  configuration  for  an  operational  amplifier  Is  shown 

ILn  Fig.  C-2  with  a  source  generator,  e  ,  and  source  resistance,  R  ,  connected. 

s  s 


Figure  C-2.  Nonlnverttng  op  aap  circuit  with  eource. 

I'he  noninverting  configuration  would  normally  be  used  with  a  high  Impedance 
6:ource  since  the  Input  Impedance  is  high.  The  circuit  model  Including  the 
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resistor  and  amplifier  noise  sources  Is  shovm  In  Fig.  C-3. 


Figure  C-3.  Inv«rtlng  op  oap  with  •ource. 

In  this  circuit,  e  can  either  be  the  noise  voltage  associated  with  R  cr  the 
s  s 

signal  voltage.  Since  all  of  the  noise  sources  (e  .  e  ,  e  ,  e  ,  1  )  are 
assumed  to  be  Independent,  consider  theli  effects  separately.  Pick  one 
source,  set  the  others  to  zero  (replace  voltage  sources  by  short  circuits, 
replace  current  sources  by  open  circuits),  and  calculate  e^  (the  output 
voltage)  for  each.  Once  this  has  been  done  for  all  of  the  noise  sources,  add 
the  squares  of  the  Individual  output  voltages  to  get  the  total  mean-square 
noise  voltage.  Calculate  e^  for  the  signal  source,  square  It,  and  divide  by 
the  total  mean-square  noise  voltage  to  get  the  output  signal-to-noise  ratio. 
This  is  also  the  equivalent  Input  slgnal-to-noise  ratio  for  the  equivalent 
Ideal  amplifier.  If  the  equivalent  Input  mean-square  noise  voltage  Is 
required,  divide  the  total  mean-cquare  noise  voltage  by  the  circuit  voltage 
gain  squared. 
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For  example,  select  e  to  analyze  and  set  all  the  other  sources  to  zero. 

5 

The  resulting  circuit  is  given  in  Fig.  C>4. 


Figure  C-4.  Circuit  of  Fig.  C-3  with  only  source  generstor  sctlve. 


Two  of  the  voltage  loop  equations  are 

i^R  +  e^  -  eg  -  0  (C-1) 

and  e-e+iR-Ae  *0  (C-2) 

S  1  2  2  1 

Eliminate  i  ,  let  A  become  very  large,  and  solve  for  Ae,  (which  equals  e  ): 

2  .  w  j 


out 


Ae 


S  2  p 


(C-3) 


where  R  is  the  parallel  combination  of  R  and  R  ,  which  equals 


Since  the  source  can  represent  either  the  signal  voltage  source  or  the 
noise  source  of  the  source  resistance,  there  are  two  results  from  this 
calculation.  One  is  the  signal  voltage  gain  of  the  amplifier: 


e  /o 
out  s 


R  /R 

2  p 


(C-4) 


The  ether  is  the  mean-square  noise  voltage  from  the  source  resistance,  which, 
at  the  output  is 


ojt-S 


4k  TR  (R  .'R  )‘ 

B  S  2  p 


(C-5) 
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and  referred  back  to  the  Input  is 


in-s 


e*  /G® 

out-S 


4k  TO 

B  S 


(C-6) 


The  other  noise-voltage  sources  are  treated  similarly. 


To  analyze  the  amplifier  noise  current,  set  all  other  sources  to  zero  to 
obtain  the  circuit  In  Fig.  C-5. 


Figure  C-S.  Circuit  of  Fig.  C-'O  with  only  aaplirier  currant  nolaa. 

Two  of  the  loop  equations  are: 

(1+1  )R  +  e  +1R  -  0  (C-7) 

n  2  1  1  n  S 

and  -IP.  -  e  +1R  -Ae  =  0  (C-8) 

n  S  i  2  2  1 

Eliminate  1^,  let  A  become  very  large,  and  solve  for  Ae^: 

e  «  Ae.  «  -  1  R,(l  +  R/B,  )  {C-9) 

out  1  n  2  S  P 

SO  the  mean-square  noise  voltage  from  the  amplifier  current  noise  Is 

,  -  [1  R  (1  +  R„/R  ))*  (C-10) 

out-l  i>  2  S  p 

at  the  output  and 

ef  .  =  1^  (R  +  R..)*  (C-11) 

In-l  n  p  S 

referred  back  to  the  input. 
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The  outpv^t  signal -to-nolse  ratio  Is  equal  to  the  voltage~aquared  signal 
output  (Eq.  C-3  squared)  dlv.'.ded  by  the  sun  of  the  voltage -squared  noise 
outputs  {twc  of  which  are  given  by  Eqs.  C-5  and  C-10),  For  all  of  the  noise 
sources,  the  output  slgnal-to~notse  ratio  for  this  circuit  Is  then 


SNR 

out 

ef  /  te^  +  4k  +  R  )  ♦  1^(R  + 

S  n  BSp  npS 


(C-12) 


This  Is  also  the  equivalent  Input  slgnal-to-nolse  ratio  and  the  denominator  Is 
the  mean-square  noise  voltage  referred  to  the  Input. 


Note:  In  general,  the  source  has  a  complex  Impedance,  not  Just  a 

resistance.  The  thermal  noise  generated  directly  by  that  Impedance  Is  only 
generated  by  the  real  part  so.  In  the  noise  terms  that  are  of  the  form  4k  TR, 

B 

only  the  real  part  of  the  source  impedance  should  be  used.  However,  In  those 
terms  that  result  from  currents  flowing  through  the  source  Impedance  (result¬ 
ing  from,  for  example,  e  or  1  ) ,  the  resulting  voltage  depends  on  the 

n  n 

magnitude  of  the  Impedance  and  so  |Z  |  should  be  used  In  place  of  R  . 

s  *  s 

Inverting  Amplifier 


Ths 

attached 


Inverting  configuration  for  an  operational  amplifier  with  a  source 
is  shown  in  Fig.  C-6. 


Figjre  C-6.  Inverting  conf  1  guretlon  for  op  anp. 
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This  configuration  would  normally  be  used  with  a  low-impedance  source  since 
the  input  Impedance  Is  low.  Figure  C-T  Illustrates  the  circuit  model 
Including  the  resistor  and  amplifier  noise  sources. 


Figure  C-7.  Equivalent  nolae  Kndet  for  Inverting  cip  acp. 

Obviously,  open-circuiting  the  output  does  not  slmplilty  the  circuit  —  It  Is 
already  open  circuited,  so,  to  proceed  with  this  techxilquc,  short  the  output 
to  ground  (and  Insert  a  small  resistor  representing  tlie  output  Impedance  of 
the  op).  The  resulting  circuit  is  shown  In  Fig.  C-EI. 


Figure  C-8.  Shor ted-output  nodel  for  Inverting  fenpllfler. 


(Notice  that  the  output  current  really  also  Includes  the  current  that  flows 

through  R  ,  but,  if  A  is  large,  this  centribution  is  negligible.) 

2 

From  this  point  on,  the  solution  procedure  is  similar  to  the  analysis  of 

the  noninverting  amplifier.  First,  compute  the  output  current,  1  ,  for  only 

0 
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the  signal  source  voltage,  e  .  Squaring  this  gives  a  quantity  proportional  to 

S 

the  signal  output  power;  this  quantity  will  be  the  numerator  of  the  slgnal-to- 
noise  expression.  Next,  compute  the  individual  components  of  noise  output 
current  by  considering  each  noise  generator  separately.  Finally,  form  the 
output  signal-to-noise  ratio  by  dividing  the  mean-square  signal  cuirent  by  the 
sum  of  the  mean-square  noise  components.  The  mean-square  signal  output 
current  for  this  circuit  is 


out 


(A/Tl^)^(e  R /R  f 

0  S  2  Bun 


(C-13) 


where  R  =  R  +  R  +  R  ,  and  the  output  signal-to-nolse  ratio  is 

•tun  S  1  2  o 


SNR 


output 


e^  /  ((R  /  R,)^ 

S  auB  2  n 


4k  TR  (1  ^  R  /R  )  +  1^  R^] 

Pa  a  2  n  a 


(C-14) 


As  in  the  case  of  the  noninverting  amplifier,  Eq.  C-14  is  also  the  input 
slgnal-to-noise  ratio  and  the  denominator  Is  the  equivalent  input  noise 
voltage.  Unfortunately,  the  process  of  referring  the  noise  to  the  input  can 
cause  some  confusion  here.  The  denominator  of  Eq.  C-14  is  referred  to  the 
positive  terminal  of  the  source  generator,  e^.  The  conventional  point  to 
which  the  equivalent  input  values  are  referred  is  the  node  between  the  source 
resistance,  R^,  and  the  resistor.  R^.  If  it  is  necessary  to  find  the  equiva¬ 
lent  voltage  referred  to  the  conventional  point,  the  analysis  should  be  done 
as  for  the  noninverting  amplifier  being  careful  to  calculate  the  circuit 
voltage  gain  from  the  proper  node  to  the  output  node.  However,  the  end  result 
for  circuit  evaluations  for  noise  is  usually  signal-to-noise  ratio;  if  that  is 
true,  thcr  it  is  easier  to  calculate  the  SNR  directly  as  shown  above. 

(Note:  The  expression  given  in  the  text  for  the  equivalent  input  noise 
voltage  of  the  inverting  amplifier  is  referred  to  the  conventional  node.  For 
the  noninverting  amplifier,  as  long  as  the  amplifier  input  impedance  is  high 
compared  to  the  source  impedance,  this  distinction  is  irrelevant.) 
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Field-Effect  Tm^l?t9rs  IJEETsJ, 


JFETs  are  sometimes  used  as  the  first-stage  preamplifier  for  sensors  with 
high  source  impedance.  The  simplified  model  [68]  for  a  JFET  shown  In  Fig.  C-9 
Is  often  adequate  for  calculating  the  noise  performance  of  a  sensor/pre¬ 
amplifier  combination. 


J' 


g  e 

■ 


Figure  C-9.  Nolae  nodal  for  JFET. 


Typical  applications  are  as  an  amplifier  or  as  a  follower  (Fig.  C-10). 


Figure  C-10.  Amplifier  and  follower  applications  for  JFETs. 


The  follower  Is  used  as  a  hlgh-lnput-impedance  buffer  to  drive  an  amplifier 
stage,  while  the  amplifier  circuit  provides  some  gain.  In  the  amplifier 
circuit,  R  sets  a  DC  bias  point  and  the  capacitor,  C  ,  effectively  removes 

B  B 

the  resistor  for  signal  frequencies. 

To  analyze  the  noise  characteristics  of  the  amplifier,  connect  a  source 
generator,  e  ,  and  its  source  resistance,  R  ,  replace  the  R.-C  combination  by 

S  5  B  B 

a  short  circuit,  short  the  output  to  ground,  and  compute  the  signal-to-noise 
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ratio  based  on  the  output  current.  The  equlvalont  circuit  for  this  analysis 
Is  given  In  Fig.  C-11. 


Figure  C-11.  Equivalent  circuit  for  noise  analyale  of  coaeton- 
eource  JFET  aaplirier. 


For  signal,  the  output  current  squared  Is 


out 


gV  (R  /R^)^ 

■  S  P  s 


(C-15) 


where  R  *  R  R /(R^  +  R..).  while,  for  noise,  the  mean-square  output  current  is 

p  S  G  S  G 


out-nolae 


g^  [e^  +  1^  R* 

m  n  n  p 


(C-16) 


+  4k„TR  +  4k„TR /(g  R„)'j 

Bp  B  D  ■  D 


For  an  amplifier,  g  R„  >  1,  and.  In  the  usual  application,  R_  «  R  ,  so  the 
last  term  in  square  brackets  can  be  neglected.  The  slgnal-to-nolse  ratio  Is 
Eq.  C-15  divided  by  Eq.  C-16.  The  equivalent  mean-square  noise  voltage 
referred  to  the  JFET  gate  Is 


e*  -  e®  +  1^  R^  +  4k  TR  (C-17) 

n  n  p  Bp 

As  discussed  previously.  If  the  source  Impedance  Is  complex,  then  Just  the 

real  part  of  R  Is  used  In  the  4k  T  term  while  the  magnitude  of  the  parallel 
P  "  2 

combination  of  the  source  Impedance  and  R  Is  used  In  the  1  term. 
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Microminiature  sensors  tend  to  be  high- impedance  devices,  hence  the 
emphasis  on  JFET  circuits.  For  a  low-impedance  sensor,  if  an  op  amp  did  not 
provide  adequate  performance,  a  bipolar  transistor  input  stage  would  be 
appropriate  [50,  54].  The  design  and  analysis  procedures  are  similar  although 
a  somewhat  more  complicated  model  than  was  used  for  the  JFET  is  usually 
required  for  the  bipolar  transistor. 
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Appendix  D.  Electronechenical  Tranaducer  Equivalents 

A  complete  analysis  of  noise  in  a  transducer  system  should  ccinsider  the 
mechanical  system  coupled  to  the  electrical  system  through  the  first  stage  of 
amplii ication.  This  Appendix  reviews  the  equivalent  circuit  repnrsentation 
for  mechanical-to-electrical  transducer  systems  of  three  types:  reciprocal, 
antireciprocal,  and  nonreciprocal.  Once  the  proper  equivalent  circuit  has 
been  drawn,  the  noise  analysis  follows  by  attaching  a  Johnson-roise  voltage 
generator  to  each  resistance  In  the  circuit. 

Both  the  reciprocal  and  the  antireciprocal  transducers  allow  transduction 
in  either  direction  (receiving  or  transmitting).  While  the  focus  here  is  on 
conversion  of  pressure  or  acceleration  to  voltage  or  current,  these  transdu¬ 
cers  will  generate  pressure  or  motion  if  a  current  or  voltage  is  applied. 

Most  conventional  transducers  fall  into  one  of  these  categories;  however, 
there  is  a  very  important  tyjse  of  transducer  —  the  electron-tunneling  sensor 
—  in  which  the  reverse  action  is  governed  by  a  different  mechanism  than  the 
forward  action.  The  reverse  action  is  so  small  that  the  transducer  can  be 
considered  to  be  practically  unidirectional  (displacement  input  produces 
current  output).  These  nonreclprocal  transducers  are  straightforward  to  model 
but,  because  a  resistance  corresponding  to  the  mechanical  damping  does  not 
translate  over  into  the  eloctrlcal  side  of  the  equivalent  circuit,  sometimes 
the  mechanical-thermal  noise  is  in:orrectly  neglected. 

The  reciprocal  and  tie  antireciprocal  transducers  can  be  represented  by 
a  "black-box"  (Fig.  D-1)  in  which  force,  F,  and  velocity,  u,  on  one  side  are 
related  to  voltage,  V,  and  current,  I,  on  the  other  side  [11]: 


Figure  D-1.  Two-port  reproientati on  for  alectroBechanical  tranaducers. 
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If  the  system  Is  linear,  the  simplest  general  relations  between  the  quantities 
are 


F  -  T  I  + 

M 

where  is  the  electrical  Impedance  with  the  iiechanical  notion  blocked 
(u  <■  0),  Z  Is  the  mechanical  Impedance  with  the  electrical  terminals  open 

i40 

(I  «  0),  T  is  the  transduction  factor  from  velocity  to  open-'circult  voltage, 

on 

and  T  Is  the  transduction  factor  from  current  to  force  when  motion  of  the 

Me 

device  Is  blocked. 

If  T  ■>  T  ,  then  the  device  Is  reciprocal.  A  transformation 
em  ae 

coefficient,  0,  can  be  defined 

^  -  T  /  Z^  (D-2) 

CIS 

and  one  way  of  drawing  the  electrical  equivalent  circuit,  which  includes  the 
mechanical“to-electrlcal  conversion,  is  shown  In  Fig.  D-2. 


Figure  D-2.  Generellzed  equivalent  circuit  for  reciprocal 
tranaducer . 

Here,  the  radiation  Impedance,  Z  is  included  separately  since  it  is  often 

rod 

negligible.  The  quantity,  Z^^,  then  refers  to  the  mechanical  Impedance  with 
the  electrical  terminals  shorted  and  In  the  absence  of  the  radiation  load. 

If  T  B  -  T  ,  then  the  transducer  is  antireciprocal.  Another 

em  me 

transformation  coefficient,  (A  ,  can  be  defined: 

H 

^  ■  T  (D-3) 

K  em 


T  u 


Z  u 

mo 


(D-l) 
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and  an  equivalent  circuit  such  as  that  in  Fig.  D-3  can  be  dravm. 


Figure  U-3.  Cenerelized  equWelent  circuit  /or  Mitl-reelproeel 
tranaducar . 

As  before  the  radiation  admittance  (Y  »  1/Z)  is  included  separately.  Notice 
that  the  mechanical  voltage-like  quantity  is  velocity  and  the  meclanical 
current-like  quantity  is  force.  This  switch  from  the  normal  impedance  analogy 
was  made  to  accomodate  the  antireciprocal  transduction  coefficients. 


For  reciprocal  transducers,  a  more  specific  circuit  (Fig.  D-4)  can  be 
drawn  that  represents  several  Important  types  [65^  69,  70]: 


r 


R 

rad 


m  /r 
rad 


R  /4> 

M 


Sp/<p 


Figure  D-4.  Typical  equivalent  circuit  for  capacitive  or 
piezoelectric  (reciprocal)  hydrophunee. 


This  circuit  can  be  used  to  analyze  simple  electrostatic  (i.e. ,  capacitive) 
and  piezoelectric  sensors.  One  of  the  assumptions  invoked  to  draw  this 
circuit  is  that  the  sensor  is  small  with  respect  to  an  acoustic  wavelength  so 
the  force  on  the  sensor  face  is  equal  to  the  acoustic  pressure,  p,  times  the 
sensor  face  area,  S.  For  hydrophones  operating  below  several  tens  of  kilo¬ 
hertz,  the  radiation  resistance  is  negligible;  however,  depending  on  the 
construction  of  the  device,  the  radiation  mass  loading  may  be  significant  (see 
Appendix  E).  For  microphones  operating  over  high  audio  and  ultrasonic 
frequencies,  the  radiation  resistance  should  be  included  [71]. 
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Capacitive  sensors 

One  of  the  simplest  reciprocal  sensors  Is  the  electrostatic  or  capacitive 
sensor  [11].  Two  parallel  plates,  each  of  area  S,  are  spaced  a  distance,  x, 
apart  and  a  polarizing  voltage  is  applied.  Vfhen  the  plate  spacing  changes 
because  of  an  applied  pressure,  the  charge  on  the  plates  changes,  which  is 
sensed  electrically.  In  the  following  discussion,  it  will  be  assumed  that  the 
polarization  circuit  has  no  if feet  on  the  noise  of  the  sensor.  (The  resistor 
—  typically  very  large  —  i.i  series  with  the  polarization  supply  is  effect¬ 
ively  in  parallel  with  so  the  electrical-thermal  noise  is  slightly 
reduced. ) 


The  capacitance,  and  the  transformation  coefficient,  are  given  by 


O' 

C 

”  C  S  /  x^ 

(D-4) 

0 

0 

»■  V  /  x^ 

0  0  0 

(D-5) 

where  c  is  the  permittivity  of  the  material  between  the  plates,  is  the 
spacing  between  the  plates  with  no  pressure  applied,  and  is  the  polariza¬ 
tion  voltage.  The  resistance,  R^,  can  be  measured  or,  if  the  loss  tangent  is 
known  for  the  dielectric,  calculated. 

The  mass  can  be  calculated  by  adding  the  mass  of  the  moveable  plate  (or 
membrane)  to  the  radiation  mass  (see  section  on  radiation  impedance).  Then 
v.he  resonance  frequency  and  Q  can  be  measured  with  the  electrical  terminals 
shorted  to  get  the  mechanical  resistance  and  the  spring  constant.  If  the 
moveable  plate  Is  a  membrane  and  the  membrane  tension  7  is  known,  then  the 
stiffness,  k,  is  8n7.  Alternately,  if  the  receiving  sensitivity  as  a  function 
of  frequency  can  be  measured,  then  the  mechanical  parameters  can  be  inferred 
from  those  measurements  combined  with  an  analysis  of  the  equivalent  circuit. 
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Piezoelectric  sensors  are  similar  to  capacitive  sensors  in  that  an 
applied  strain  produces  a  charge  redistribution.  There  are,  however,  many 
more  configurations  in  which  pie:':oelectric  material  can  be  used,  the  choice 
depending  on  the  application.  Here,  three  common  configurations  will  be 
considered:  a  hydrostatic-mode  sensor,  a  cylinder,  and  a  bender-disk.  The 
values  in  the  equivalent  circuit  will,  in  general,  be  different  for  each 
configuration  but  the  form  of  the  circuit  will  be  the  same  or  very  close  to 
that  used  for  the  capacitive  sensor. 

The  hydrostatic-mode  sensor  is  simply  a  rectangular  block  of  plezoceramlc 
with  electrodes  on  one  pair  of  opposing  faces  [65].  Pressure  acts  uniformly 
on  all  sides,  the  mechanical  response  is  stiffness-dominated,  and  the  device 
can  tolerate  very  large  pressure  loads.  Unfortunately,  this  configuration  has 
rather  low  sensitivity. 

The  same  equivalent  circuit  can  be  used  as  was  used  for  the  capacitance 
sensor  (and  the  mechanical  mass  can  be  neglected).  The  capacitance,  C^,  and 
the  transformation  coefficient,  are  given  by 


C^3  S  /  t 

(D-6) 

^0^3 

(D-7) 

where  is  the  dielectric  constant  for  zero  strain,  t  is  the  block  thickness 

33 

from  electrode  to  electrode,  S  is  the  surface  area  of  one  electrode,  and  h^^ 
is  the  rate  of  change  of  the  electric  field  with  applied  strain  for  constant 
electric  flux  density.  (The  33  subscript  refers  to  the  relationship  between 
the  direction  of  the  field  from  electrode  to  electrode  and  the  direction  in 
which  the  material  was  polarized:  For  the  33-mode,  these  directions  are 
the  same.  These  parameters  can  be  found  in  numerous  references  on 
piezoelectric  materials,  for  example  (65,  72].) 
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The  mechanical  stiffness,  k  ,  is 

k  =  S  /  [t(s,,  +  2s„)I  (D-8) 

n  11  12 

where  the  s^^  and  s^^  are  piezoceramic  elastic  constants.  Since  the  transdu¬ 
cer  is  very  stiff,  the  radiation  mass  loading  would  normally  be  negligible; 
since  this  configure t,ion  is  normally  only  used  as  a  hydrophone,  the  radiation 
resistance  loading  can  be  neglected  as  well.  Under  these  conditions,  the 
open-circuit  receiving  sensitivity  is 

where  the  g’s  are  the  rates  of  change  of  electric  field  with  respect  to 
changes  in  applied  stress  for  constant  electric  flux  density.  (Usually,  g^^ 
is  about  half  of  g^^  and  of  the  opposite  sign  so  the  “hydrostatic"  sensitivity 
is  low. )  If  the  loss  tangent  of  the  piezoceramic  material  is  known,  then 
can  be  calculated. 

Higher  sensitivities  can  be  obtained  by  isolating  one  or  more  faces  of 
the  material  from  the  acoustic  pressure.  This  is  often  accomplished  by  using 
an  air  cavity.  One  such  configuration  is  a  cylinder  of  piezoceramic  material 
with  end  caps  and  air  Inside,  for  a  cylinder  of  length,  h,  radius,  a,  and 
thickness,  t,  the  mechanical  mass.  m.  is  p2iiaht  where  p  is  the  density  of  the 
ceramic.  The  applied  force,  F,  is  equal  to  the  acoustic  pressure  times  the 
surface  area,  2nah. 

If  the  cylinder  is  polarized  parallel  to  the  cylinder  axis  and  the 
electrodes  are  on  the  inner  and  outer  surfaces  of  the  cylinder,  thri  the 
device  is  being  operated  in  the  3-1  mode  [65,  73].  The  capacitance,  C^,  in 
this  case  is 

C  =  e®2nah/t  (D-10) 

0  31 

and  the  mechanical  stiffness  is 

k  =  2iiht/aB^  (D-ll) 

m  n 
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The  transformation  coefficient  Is 

2  It  h  /  (D-12) 

and  the  mechanical  resistance  is  best  determined  by  measuring  the  Q  with  the 
electrical  terminals  shorted. 

An  even  simpler  configuration  that  it.  air-  is  the  bender  disk.  A 

thin  piszoceramlc  disk  Is  cemented  to  a  brass  disk  that  is,  in  turn  edge 
mounted  over  an  air-filled  cavity  [73] .  Most  often,  two  such  disks  arr 
mounted  back-to-back  with  an  edge  ring  to  separate  them  and  ’.o  ^orlrt  the 
cavity.  If  the  brass  disk  (radius,  a,  and  thickness,  t)  dominates  ^.udss 
and  stiffness,  then  the  mechanical  stiffness  Is 

k  =  37ift^/2a^  (D  '13) 

n 

where  Y  is  the  Young’s  modulus  of  the  brass  disk,  and  the  mechanical  mass  is 

m  =  2Tta^tp/3  (D-14) 

where  p  is  the  density  of  the  disk.  The  effective  water  mass  should  be  added 
to  get  the  total  mass.  This  additional  mass  is 

m  =  0.4  n  a^  p  (D-15) 

w  Hater 

The  capacitance,  C^,  is 

C  =  S  /  t  (D-16) 

0  33 

Since  the  disk  does  not  move  as  a  rigid  piston,  Eq.  D~14  is  less  than  the 
total  mass  of  the  disk.  This  also  complicates  the  calculation  of  ^  so  that 
quantity  should  be  obtained  Irom  the  measured  pressure  response  of  the  device. 
The  losses  can  be  determined  from  the  measured  Q. 

Moving-coil  sensors 

The  moving-coil  sensor  is  the  archetypical  antireciprocal  transducer. 
Because  the  output  is  proportional  to  the  velocity  of  the  proof  mass  rather 
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than  Its  displacement,  it  Is  more  convenient  to  use  the  mobility  analog 

(velocity  i - »  voltage;  force  < - >  current)  lor  the  mechanical  parts  of  the 

circuit.  A  generic  equivalent  cl^.'cuit  for  those  BOvlr*g-ccll  devices  (for 
example,  [22]}  is  shovm  in  Fig.  D-5. 


rigure  D-5.  Equivalent  circuit  for  a  ■nvtng-coll  eenBor. 

The  mass-spring  oscillator  is  represented  by  m,  R  ,  and  k  ;  the  electrical 

n  B 

properties  of  the  coil  by  H  and  L  ;  and  the  mass  of  the  sensor  case  by  m  . 

c  c  c 

The  shunt  resistor  is  used  to  adjust  the  damping  of  the  sensor  and  is  usually 
ter-  il  i  the  radiation  load  is  not  negligible,  it  can  be  added  in 

parai  x  with  the  velocity  ge;  ~'^ator  as  in  Fig.  D-6. 


Figure  D«<’6.  Addition  of  radiation  load  for  antl-*raclprocal  aenaor. 

If  the  sensor  case  moves  at  the  same  speed  as  the  medium  in  which  it  is 
placed,  then  the  velocity  source  is  appropriate.  This  would  be  true  for  a 
small  sensor  attached  to  a  massive  object  (a  geophone  in  the  ground  or  an 
accelerometer  on  a  piece  of  machinery)  or  for  a  sensor  that  is  neutrally 
buoyant  and  not  near  resonance  (an  accelerometer  being  used  as  a  hydrophone). 
The  more  general  case  of  an  Immersed  sensor  is  better  treated  by  using  a  force 
generator  proportional  to  the  net  force  resulting  from  the  pressure  integrated 
over  the  sensor  case  (see  Appendix  E). 
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The  transformation  coefficient  is 

»  B  1  (D-17) 

2 

where  B  is  the  magnetic  flux  density  (webers/m  )  and  1  is  the  length  of  wire 
in  the  coil.  Normally,  B  is  rather  hard  to  determine;  fortunately,  can  be 
measured  [22]: 

(1)  Fix  the  case  to  a  nonmoving  surface  (u  «  0). 

(2)  Remove  the  shunt  resistor,  R  _ 

SHUKT 

(3)  Measure  the  open-circuit  damping,  b^.  The 
mechanical  resistance,  R  ,  can  be  determined  from 

fli 

(4)  Reconnect  the  shunt  resistor. 

(5)  Measure  the  damping,  b^. 

(6)  If  the  coil  inductance,  I,^,  can  be  neglected,  then 
the  transformation  coefficient  is 

.  2  »  (R^  *  R^^)  (b^  -  (D-18) 

Electron- tunneling  sensors 

For  most  transducers,  the  magnitude  of  the  output  voltage  (or  current)  is 
proportional  to  the  amount  of  sensing  material.  The  capacitive  sensor’s 
output  is  proportional  to  the  size  of  the  plates,  the  output  of  a  piezoceramic 
sensor  is  proportional  to  the  amount  of  ceramic,  and  the  output  of  a  moving- 
coil  sensor  is  proportional  to  the  amount  of  wire  i  the  ?ift.iv:?tic  flux 
density.  Miniaturization  necessarily  reduces  the  senb'tivity  of  these 
transducers. 

An  alternate  means  of  sensing  displacement  is  by  measuring  the  electron 
tunneling  current  that  flows  between  two  conductors.  If  two  conductors  are 
very  close  together  (on  the  order  of  a  nanometer),  then  there  is  a  wave- 
mechanical  transfer  of  electrons  from  one  to  the  other  [58].  This  transfer  is 
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exponentially  sensitive  to  the  separation  and  so  provides  the  basis  for  an 
extremely  sensitive  displacement  measurement.  The  surface  area  through  which 
the  tunneling  takes  place  need  only  be  a  few  atoms  in  size;  a  larger  area  does 
not  increase  the  sensitivity  to  displacement.  In  this  sense,  the  electron** 
tunneling  sensor  is  an  ideal  candidate  for  microminiaturization  [9}.  The 
mechanical  parts  of  the  accelerometer  or  pressure  configuration  are  still 
subject  to  thermal-fluctuation  motion  though. 

One  interesting  aspect  of  the  electron-tunneling  sensor  is  that  it  is 
effectively  unidirectional.  A  displacement  change  produces  a  change  in  the 
tunneling  current  but  an  externally  forced  change  in  the  current  produces  very 
little  sensor  motion.  The  motion  produced  results  from  electrostatic  forces 
on  the  equivalent  capacitance  of  the  tunneling  gap  (and  this  capacitance  is 
quite  small)  and  fr^m  momentum  transfer  from  the  tunneling  electrons  [741. 
Since  these  mechanisms  are  not  the  same  as  the  mechanism  for  the  displacement- 
to-voltage  transduction,  this  device  is  nonreciprocal.  A  convenient 
equivalent  circuit  for  this  device  is  shown  in  Fig.  D-7. 
m 


Figure  D-7.  Equivalent  circuit  for  electron-tunneling  aeneor. 


Notice  that  the  mechanical  resistance  does  not  appear  in  the  output  branches: 
The  thermal  noise  from  this  resistance  does  enter  the  output  circuit  through 
f  .  This  would  present  a  serious  problem  if  the  classical  definition  of  noise 

n 

figure  were  forced  on  this  device  since  only  R  would  be  considered  (and  R 

9  9 

is  noiseless!}.  An  advantage  of  the  nonreciprocal  nature  is  that  the  ampli¬ 
fier  current  noise  does  not  feed  back  into  the  mechanical  circuit  so  the 
amplifier-induced  noise  is  not  as  large  as  it  could  be  for  a  reciprocal 
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transducer.  This  advantage  is  Inconsequential  thoxigh,  If  the  device  Is 
dominated  by  thermal  noise  in  the  mechanical  circuit. 

In  the  equivalent  circuit  given  above. 

0  -  I  /  w*  m  (D>19) 

0  0  0 

where  I  is  the  static  tunneling  current  (a  function  of  the  bias  voltage,  V^, 
across  the  gap)  and  is  a  characteristic  displacement  associated  with  the 
tunneling  [58]: 

X  =  1  /  a  (D-20) 

0 

where  a  «  1.025  x  lo'°  m"^  and  ♦  is  the  work  function  of  the  surfaces 

between  which  the  tunneling  is  taking  place  (about  half  an  electron  volt  for 
gold-to-gold  tunneling).  The  distance,  x^.  is  usually  close  to  an  angstrom. 
Because  tunneling  Involves  random  emission  of  electrons  across  the  gap,  a  shot 
noise  component,  ,  must  also  be  Included  in  the  analysis  along 

with  the  equivalent  gap  resistance,  (Be  aware  that  tunneling  is 

not  an  equilibrium  process:  there  is  no  mechanism  to  force  the  tunneling  waves 
into  equilibrium  as  they  cross  the  gap  so  there  is  no  thermal  noise 
associated  with  R  . )  Since  this  device  is  a  current  generator,  the  first 

9 

stage  of  preamplification  is  usually  a  high-gain  current-to-voltage  convertor 
(an  Inverting  op  amp  with  R^  »  0). 
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Appendix  E.  Radiation  Loading  of  Transducers 

The  two  basic  types  of  acoustic  sensors  —  pressure  sensors  and  accelero¬ 
meters  —  differ  significantly  In  their  interaction  with  the  medium.  Most 
practical  receiving  transducers  are  small  with  respect  to  a  wavelength  in  the 
medium  over  most  of  their  operating  range.  If  this  is  true,  then  the  presence 
of  the  transducer  has  little  effect  on  the  acoustic  waves  that  pass  by.  For  a 
pressure  sensor,  this  means  that  the  relevant  pressure  is  equal  to  the 
undisturbed  acoustic  pressure  In  the  surrounding  medium. 

An  accelerometer  responds  to  case  motion  rather  than  pressure  on  the  case 
so  the  appropriate  quantity  Is  the  force  exerted  on  the  case  by  the  acoustic 
wave.  This  force  Is  equal  to  the  Integral  of  the  pressure  times  the  vector 
element  of  surface  on  the  transducer  case  over  the  entire  case.  For  a  case  of 
volume,  V,  with  dimensions  much  smaller  than  a  wavelength,  the  root-mean- 
square  value  of  the  force  resulting  from  a  traveling  plane  wave  Is 

F  =  1/  u)  p  /  c  (E-1) 

where  p  Is  the  root-mean-square  acoustic  pressure  and  c  Is  the  sound  speed  In 
the  medium. 

If  the  transducer  Is  well  represented  by  a  rigid  body  with  the  same 
density  as  the  fluid  In  which  It  is  Immersed  (that  is.  If  the  transducer  Is 
neutrally  buoyant  as  Is  common  for  underwater  sensors),  then  It  will  move  as 
If  It  were  a  part  of  the  fluid.  By  using  the  above  expression  for  force  with 
the  proper  mechanical  equivalent  circuit,  the  proper  case  motion  can  be 
computed  even  If  the  density  of  the  transducer  is  not  equal  to  that  of  the 
surrounding  fluid  or  if  the  transducer  does  not  act  like  a  rigid  body.  (At 
resonance,  If  the  proof  mass  is  of  the  same  order  as  the  case  mass,  the  motion 
can  be  quite  different  even  if  the  device  is  neutrally  buoyant. ) 

Another  difference  between  the  pressure  sensor  and  the  accelerometer  is 
in  the  radiation  impedance.  The  pressure  sensor  acts  as  a  monopole  source  as 
far  as  reradlatlon  is  concerned,  while  the  accelerometer  acts  as  a  dipole 
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source  for  reradlatlon.  A  reasonable  model  of  a  small  pressure  transducer 
(for  the  purposes  of  radiation  Impedance  calculation)  Is  as  an  unbaffled 
piston.  In  this  case,  the  (mechanical)  radiation  Impedance  is  [11] 

Z  “  p  c  n  a^  fdca)*  /  4  ♦  1  0.6  kal  (E-2) 

MECH  I  I 

whore,  here,  k  is  the  acoustic  wave  number  (w/c)  and  a  is  the  radius  of  the 
piston.  (As  in  the  body  of  this  paper,  the  convention  e  Is  used  for 
convenience  In  drawing  electrical  equivalent  circuits. ) 

A  simple  model  for  the  reradlatlon  of  an  accelerometer  Is  that  of  a 
sphere  oscillating  along  some  axis.  The  radiation  impedance  corresponding  to 
this  motion  is  [30] 

Z  »  p  c  n  a^  [(ka)*  /  3  +  1  2  ka  /  31  (E-3) 

MECH  L  J 

In  comparing  Eqs.  £-2  and  £-3,  notice  that  the  Imaginary  parts,  representing 
the  mass  loading,  are  roughly  equal;  however,  the  real  parts,  representing  the 
real  radiation  and  contributing  to  the  thermal  noise,  are  significantly 
different.  The  dipole  radiation  is  proportional  to  (ka)*  which,  for  small  ka. 
Is  much  smaller  than  the  monopole  dependence  of  (ka)  . 

If  the  accelerometer  Is  neutrally  buoyant,  the  noise  associated  with 
reradlatlon  resistance  can  be  calculated  as  an  equivalent  pressure  as  is  done 
In  the  text.  This  noise  component  can  then  be  treated  as  an  addition  to  the 
ambient  noise.  If  the  accelerometer  Is  not  neutrally  buoyant  (an  accelero¬ 
meter  In  air  would  probably  not  be),  then  Eq.  E-3  should  be  used. 
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